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SUMMARY 


The  measurement  of  ln-situ  stress  la  important  to  many  OKA  ground 
notion  and  atructuraa  programs.  Recent  result a  from  field  teats  have 
shown  that  long  duration,  reproducible  records  can  be  obtained  with 
armored  gages  over  a  very  large  range  of  stress  and  In  severe  mechanical 
and  M  environments.  However,  reduction  of  the  records  to  stress  data 
requires  that  the  relationship  between  free-field  stresa  and  stress  In 
the  gage  and  Its  Immediate  surroundings  (borehole)  be  established  for 
the  variety  of  geologic  materials,  borehole  fillers,  and  stress  ranges 
encountered  In  DMA  field  tests.  Although  such  relationships  have  been 
Investigated  both  analytically  and  experimentally  for  elastic  response 
of  materials,  the  elastic-plastic  and  plastic  response  reglaea  have 
received  little,  If  any,  attention.  These  regimes  are  of  particular 
Importance  In  DMA  programs  because  the  yield  behavior  of  in-sltu 
geologic  materials  subjected  to  particular  strain  paths  Is  largely 
unknown  and  Is  difficult  to  determine  in  the  laboratory.  Recent 
advances  In  finite  difference  and  finite  element  calculatlonal  simula¬ 
tions  have  provided  new  tools  for  Investigating  the  required  atrese 
relationships. 

SRI  Is  conducting  a  theoretical,  calculatlonal,  and  experimental 
program  to  determine  the  relationships  between  Inclusion  (gage,  bore¬ 
hole)  stress  and  free-field  (medium)  stress  for  materials,  geometries, 
and  stress  ranges  encountered  In  ln-sltu  measurement  programs  and  to 
provide  criteria  for  gage  design  and  installation.  The  program  consists 
of  two  complementary  efforts:  (1)  finite  element  calculations  with 
supporting  analytical  computations  and  static  trlaxial  compression 
experiments,  and  (2)  finite  difference  calculations  with  supporting 
shock  wave  experiments  to  examine  dynamic  effects.  In  the  first  effort, 
we  ere  addressing  conditions  encountered  when  stresses  of  wavelengths 
that  are  large  compared  with  Inclusion  dimensions  are  to  be  measured. 


and  in  th«  second,  we  are  addressing  stresses  of  wavelength  comparable 
to  the  Inclusion  dlwenslon. 

We  have  performed  preliminary  biaxial  plane  strain,  finite  element 
calculatlona  to  determine  the  stress  distribution  normal  to  gage  planes 
in  a  borehole,  in  tuff,  filled  with  a  matching  grout  and  loaded 
proportionally*  We  have  investigated  the  effect  of  material  property 
mismatch,  grout-tuff  bonding  (including  separation),  and  unloading 
(residual  stress)  on  stress  distribution  and  ratio  of  borehole  stress  to 
free-fleld  stress.  The  preliminary  calculations  were  restricted  to 
elastic  free-field  loading,  but  they  allowed  plastic  behavior  in  the 
inclusion  and  in  the  medium  next  to  the  inclusion.  These  results 
indicated  that: 

(1)  The  dletrlbtttion  of  atreaaea  normal  to  the  radial  and 
tangential  atreaa  measurement  planes  in  the  borehole  is 
nonuniform.  The  nonuniformity  is  a  strong  function  of 
applied  load  and  interface  bonding. 

(2)  Over  the  central  portion  of  the  radial  meaauring  plane 
(50Z  of  borehole  diameter),  the  distribution  is  uniform 
within  <102  for  bonded,  frictional,  and  unbonded  inter- 
facea,  but  the  ratio  of  borehole  streas  to  free-field 
atreaa  is  closer  to  unity  for  bonded  and  frictional 
Interfaces. 

(3)  The  stress  distributions  on  the  tangential  measuring 
plans  warm  more  uniform  for  the  bonded  Interface; 
however,  the  stress  ratio  depended  on  the  free-field 
stresa  amplitude. 

(4)  Separation  at  the  interface  cauaad  nonuniform  stress 
distributions  similar  to  those  caused  by  slip. 

(5)  kesldual  stress  remained  in  the  borehole  after  unloading 
of  the  free-field  stress.  The  magnitude  and  dlstrlbu- 
tlona  were  dependent  on  the  boundary  conditions,  match¬ 
ing  of  properties,  and  orientation  of  tha  measurement 
plane.  The  least  residual  stress  and  most  uniform 
distribution  occurred  for  the  radial  maaauramant  plane 
and  frictional  boundary.  The  bonded  boundary  produced 
large  residual  atresa  and  larga  variations  in  distribu¬ 
tion  if  die  properties  were  not  well  matched;  therefore, 
the  bonded  boundary  was  best  for  loading  but  poor  upon 
unloading  if  the  properties  are  not  well  matched. 


W«  hay*  also  performed  elssClc  analyses  of  stresses  in  elliptic  and 
circular  inclusions  under  plane  strain  loading  to  validate  the  finite 
element  code  and  facilitate  parametric  investigations  in  the  elastic 
regiae .  The  analyses  addressed  a  bonded  general  ellipse,  a  bonded 
slender  ellipse,  a  bonded  and  unbonded  circle,  circular  interface 
separation,  a  rigid  line  inclusion,  and  residual  stress  in  a  bonded 
circle.  Of  particular  Interest  is  the  result  of  the  rigid  line 
inclusion,  which  showed  that  the  line  (gage)  should  not  be  bonded  to  the 
medium. 

A  finite  eleaent  analysis  was  asde  of  a  typical  flatpack  stress 
gage  embedded  in  grout.  Results  showed  that  the  thin  sensing  eleaent 
was  subjected  to  the  free-fleld  stress  if  the  eleaent  occupied  less  than 
60X  of  the  central  aidplane  of  the  plastic  within  the  steel  araor. 

We  designed  and  constructed  a  static  loader  for  subjecting  borehole 
inclusions  in  10  s  10  i  22  inch  samples  to  triaxial  loads  with  a  aaximua 
stress  difference  of  ~l/2  kbar  at  a  hydrostatic  pressure  of  ~1  kbar. 
Sensors  for  these  experlaents  have  been  examined  for  suitability 
(minimizing  inclusion  effects  and  monitoring  Interface  behavior). 

We  also  performed  a  numerical  simulation  of  the  dynamic  (shock) 
loading  of  an  ytterbium  foil  In  FMMA  (experlaents  conducted  by  Y.  M. 
Gupta  at  Washington  State  University).  The  simulation  generated  the 
stress  and  strain  states  in  the  foil  during  the  load/unload  process. 

The  results  indicated  that: 

(1)  The  aaterial  models  used  for  FMMA  and  ytterbium  and  the 
plesoreslstance  formulation  of  Gupta  led  to  a  resistance 
change  versus  axial  stress  profile  in  fair  agreeaent 
with  Gupta’s  experiaent  that  measured  the  rising  and 
peak  stress  portions  of  the  loading. 

(2)  Separation  between  foil  and  PMMA  on  unloading  resulted 
from  the  geometry  of  the  foil  and  aaterial  model 
(yielding)  of  PMMA.  A  change  in  foil  geometry  on 
separation  produced  a  residual  resistance  in  the  foil. 

(3)  The  calculational  method  of  determining  the  scalar 
resistance  change  in  the  deformed  foil  should  be 
examined  further. 
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(4)  Additional  comparisons  should  be  made  between  calculated 
and  measured  response  at  other  stresses  because  adjust¬ 
ments  in  plastic  work  terms  in  the  piezoresistance 
formulation  were  required  to  obtain  quantitative 
agreement  between  measured  and  calculated  responses. 

Finally,  we  examined  the  feasibility  of  using  the  4-inch-diameter 
gas  gun  and  reduced-scale  flatpack  gages  to  investigate  the  relationship 
between  a  known  free-fleld  stress  history  and  gage  stress  history.  The 
latter  was  examined  by  performing  one-dimensional  calculatlonal  simula¬ 
tions  of  a  gage  in  FNMA.  Results  Indicated  that  gage-medium  equilibrium 
could  be  obtained  during  loading  durations  achievable  in  the  4-inch- 
dlameter  gun  only  if  the  steel  armor  of  the  gage  were  reduced  to  very 
thin  foils.  Because  this  reduction  would  alter  gage  response,  it  was 
concluded  that  longer  duration  pulses  than  are  achievable  In  the  4-lnch- 
diameter  gun  would  be  required.  This  conclusion  evoked  the  Idea  of  a 
2-ft-dlameter  flyer  plate  for  the  dynamic  experiments.  The  design  of 
this  facility  is  presented  In  this  report. 


PREFACE 
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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement . 


To  Convert  Proa 

To 

Multiply  By 
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1.000  000  X  E  -10 

atmosphere  (normal ) 

kilo  pascal  (kPa) 

1.013  2S  X  E  *2 

bar 

kilo  pascal  (kPa) 

1.000  000  X  E  *2 

barn 

aster2  (n2) 

1.000  000  X  E  -2B 

British  tharnal  unit  (thermochaaical) 

joule  (J) 

1.0S4  3S0  IE*! 
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giga  becquerel  (GBq)t 

3.700  000  X  E  *1 

degree  Celsius! 

degree  kalvin  (X) 

tK  ■  t*c  *  273.  IS 

degree  (angle) 

radian  (rad) 

1.74S  329  X  E  -2 

degree  Fahrenheit 

degree  kelvin  (X) 

tK  •  (t*F  ♦  459.67)/l .8 

electron  volt! 

joule  (J) 

1.602  19  X  E  -19 

erg! 

joule  (J) 

1.000  000  X  E  -7 

erg/second 

watt  («) 

1.000  000  X  E  -7 

foot 

mater  (a) 

3.048  000  X  E  -1 

foot -pound-force 

joule  (J) 

1.3SS  818 

gallon  (U.S.  liquid) 

meter2  (n2) 

3.78S  412  X  E  -3 

inch 

mete.-  (a) 

2.S40  000  X  E  -2 

jerk 

joule  (J) 

1.000  'OC  X  f  .9 

joule/kilograa  (J/kg)  (radiation 
dose  absorbed)! 

gray  (Gy)* 

1.000  000 

kilotonsi 

tarajoules 

4.183 

kip  (1000  lbf) 

newton  (N) 

4.44*  22 2  IE  »J 

kip/inch2  (ksi) 

kilo  pascal  (kPa) 

6.894  7S7  IE*! 

ktap 

newton-second/m2  (N-s/n2) 

1.000  000  I  E  *1 

micron 

meter  (a) 

1.000  000  X  B  -6 

ail 

meter  (m) 

2.S40  000  I  E  •! 

nile  (International) 

meter  (m) 

1 .60S  344  X  E  *3 

ounce 

kilogram  (kg) 

2.834  9S2  X  E  -2 

pound-force  (lbf  avoirdupois) 

newton  (N) 

4.448  222 

pound-force  inch 

nevton-aeter  (N*n) 

1.129  848  X  E  -1 

pound-force/ inch 

newton/meter  (N/m) 

1.7S1  268  X  E  *2 

pound -force/ foot ^ 

kilo  pascal  (kPa) 

4.788  026  X  E  -2 

pound-force/inch2  (psi) 

kilo  pascal  (kPa) 

6.894  7S7 

pound-mass  (lba  avoirdupois) 

kilogram  (kg) 

4.S3S  924  X  E  -1 

pound-mass-foot2  (moment  of  inertia) 

kilogram-meter2  (kg*m2) 

4.214  011  X  E  -2 

pound-mass/foot2 

kilogram-meter2  (kg/m2) 

1.601  846  X  E  *1 

rad  (radiation  dose  absorbed)! 

gray  (Gy)* 

1.000  000  X  E  -2 

roentgen! 

coulonb/kilograa  (C/kg) 

2.S79  760  X  E  -4 

shake 

second  (s) 

1.000  000  X  E  -8 

slug 

kilogram  (kg) 

l .459  390  X  E  ♦! 

torr  (an  Hg,  0*  C) 

kilo  pascal  (kPa) 

1.333  22  X  E  -1 

•The  gray  (Gy)  is  tha  accepted  SI  unit  equivalent  to  the  energy  i apart ed  by  ioniting  radiation  to  a  mass  of 
energy  corresponding  to  one  joule/kilograa. 

tThe  becquerel  (8q)  is  the  SI  unit  of  radioactivity;  1  8q  *  1  event/s. 

(Temperature  aay  be  reported  in  degree  Celsius  as  well  as  degree  kelvin. 

{These  units  should  not  be  converted  in  DMA  technical  reports;  however,  a  parenthetical  conversion  is 
permitted  at  the  author's  discretion. 
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SECTION  1 


INTRODUCTION  AND  PROGRAM 


1.1  INTRODUCTION 

An  important  sensor  In  Department  of  Defense  research  and  development 
field  experimenta  is  the  piexorea latent  stress  gage.  The  gage  is  uaed  to 
meaeure  the  stress  waves  in  rocks  and  soils  resulting  from  underground 
nuclear  explosions  and  from  high  explosive  systems  that  simulate  shock  and 
blast  effects  from  surface  and  aboveground  nuclear  explosions. 

One  of  the  difficulties  in  measuring  stresses  in-situ  is  that  the 
gage  package  and  its  immediate  surroundings  perturb  the  free  field.  The 
gage  package  may  be  regarded  as  an  inclusion  in  a  geologic  medium.  In 
many  cases,  the  installation  procedure  Introduces  an  additional  inclusion. 
For  example,  gages  are  installed  in  Nevada  Test  Site  (NTS)  tuff'  by 
drilling  from  a  tunnel.  Inserting  the  gage,  and  packing  the  borehole 
with  a  rock-matching  grout,  as  Illustrated  in  Figure  1.  Because  it  is 
impossible  to  match  all  the  properties  of  the  tuff  in  the  neighborhood 
of  the  gage  and  to  ensure  perfect  bonding  with  the  tuff,  the  grout  forms 
a  cylindrical  inclusion.  Thus,  the  gage  is  in  general  not  subjected  to 
the  same  stresses  as  the  free-fleld  stresses,  and  we  must  determine  the 
relationship  between  the  grout  inclusion  and  geologic  medium  stress 
fields. 

Our  coordinated  theoretical  and  experimental  program  to  investigate 
the  Inclusion  effects  of  gages  and  their  Installation  is  described  below. 
This  program  is  regarded  as  one  phase  of  the  overall  program  on  in-situ 
stress  measurement.1’  The  overall  objective  is  to  guide  the  development 
of  gage  design  and  installation  procedures  to  ensure  that  the  gage 
records  are  useful  approximations  of  the  free-field  stresses. 

The  approach  is  first  to  perform  calculations  on  current  stress 
gage  design  concepts  and  on  installation  configurations  to  identify  the 
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Influential  peranetere  with  a  view  to  optimising  the  stress  measurement 
technique  or  establishing  Its  limitations*  Second,  simple  static  and 
dynamic  experiments  are  performed  to  validate  predictions  and  to  assist 
technique  development* 

1.2  PS OGHAM 

The  Initial  phase  of  our  overall  program  on  in-aitu  stress  measure¬ 
ments  consists  of  the  following:  (1)  finite  element  code  calculations 
on  stress  gage  Inclusions  end  borehole  inclusions,  (2)  classical  anal¬ 
yses  of  these  Inclusions  at  low  stresses,  (3)  finite  difference  code 
calculations  on  the  response  of  an  ytterbium  sensing  element  to  a  planar 
shock  wave,  (4)  static  experiments  on  gages  and  boreholes,  and  (5)  e 
study  of  experimental  designs  for  providing  dynamic  loading  of  stress 
gages  and  boreholes*  A  description  of  the  theoretical  and  experimental 
portions  of  the  program  follows. 

Theoretical  Investigations 

Because  All  stress  gage  packages  fora  an  inclusion  In  the 
surrounding  medium,  it  can  be  expected  that  the  relationship  between 
gage  stress  and  free-fleld  stress  will  be  influenced  by  the  following: 

•  Bonding  of  the  gage  package  to  the  medium. 

•  Material  property  differences  between  the  gage  package  and 
medium. 

•  Package /medium  separation. 

•  Residual  stress  in  the  gage. 

•  Gage  geometry. 

Similarly,  if  the  gage  Is  placed  in  a  borehole,  the  borehole  is  an 
inclusion  in  the  surrounding  medium  and  the  stress  in  the  borehole 
depends  on  the  following: 

•  Bonding  of  the  borehole  filler  to  the  medium 

•  Plller/medium  properties  mismatch 

•  Plller/medium  separation  conditions 

•  Orientation  of  the  gage  (stress  component  measurement) 
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•  Filler  residual  stress 

•  Borehole  geometry. 

The  sensitivity  of  the  inclusion  stresses  to  the  above  effects  are 
investigated  by  finite  element  code  calculation  (Sections  2  and  4).  For 
■any  field  applications,  the  stress  gradient  is  saall  compared  with  the 
inclusion  dimensions,  so  static  solutions  are  appropriate  for  these 
cases. 

Solutions  have  been  generated  by  classical  analysis  for  cases  where 
the  Inclusion  and  medium  remain  elastic,  except  for  residual  stress 
cases  where  the  filler  (grout)  has  been  allowed  to  undergo  plastic 
deformation  (Section  3).  These  closed- form  solutions  are  useful  for 
parametric  studies  and  the  validation  of  finite  element  codes  in  the  low 
stress  regime. 

The  other  main  area  of  theoretical  investigation  in  this  program  is 
the  response  of  an  ytterbium  foil  sensing  element  to  shock  waves 
(Section  5).  This  work  has  been  performed  by  using  a  finite  difference 
wave  code,  and  it  forms  support  for  the  basic  experimental  and  theoretical 
work  on  piezoresistance  being  conducted  for  DMA  by  T.  M.  Gupta.  The 
calculations  correspond  to  gas  gun  experiments  in  which  ytterbium  foil 
embedded  in  PMMA  was  subjected  to  a  planar  shock  wave  generated  by  an 
impacting  flyer  plate. 


Experimental  Investigations 

The  experimental  program  comprises  two  parts:  (1)  stress  gage 
package  design  and  performance,  and  (2)  borehole  inclusion  effects.  For 
the  stress  gage  the  work  consists  of  the  following: 

•  Gage  design  and  fabrication 

•  Loading  and  unloading  of  gage  In  fluid,  sand,  grout 
(static),  and  PMMA  (dynamic). 

The  loading  and  unloading  experiments  are  designed  to  provide  deter¬ 
minate  or  known  stress  environments  for  comparison  with  the  stress  gage 
measurements,  including  possible  residual  stresses  Induced  by  the  gage 
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Inclusion*  Appropriate  results  will  be  correlated  with  theoretical 
predictions* 
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For  the  filled  borehole  Inclusion  an  experiment  has  been  designed 
end  the  apparatus  constructed,  as  described  In  Section  6.  In  the 
experiment,  e  general  trlaxlal  principal  stress  field  can  be  established 
la  a  rectangular  block  of  geologic  material  or  rock  simulant*  A  filled 
borehole  can  be  provided  along  the  major  axis  of  the  block  so  that 
stress-strain  measurements  can  be  compared  in  three  basic  cases:  (1)  no 
borehole,  (2)  borehole  with  bonded  grout  fill,  and  (3)  borehole  with 
unbonded  grout  fill*  The  experiment  Is  also  suitable  for  parametric 
studies  (strengths,  loading  ratios,  separation,  geometry,  and  so  forth). 
Suitable  sensors  are  examined. 

Dynamic  (gas  gun)  loading  of  flatpack  gages  was  investigated  by 
numerical  simulation  to  determine  the  feasibility  of  reducing  gage 
dimensions  appropriate  to  the  dimensional  limitations  of  the  4-lnch- 
diameter  gas  gun.  Because. this  Investigation  showed  that  the  gas  gun 
cannot  provide  pulses  of  adequate  duration  for  flatpack  gages,  the 
design  of  a  larger  flyer  plate  facility  was  added  to  the  program 
(Section  6). 


SECTION  2 


STATIC  CALCULATIONS  FOR  FILLED  BOREHOLE  INCLUSIONS 

2.1  INTRODUCTION 

Finite  element  calculations  have  been  made  for  the  filled  borehole 
Inclusion  to  let  us  coapare  the  stresses  acting  normal  to  the  gage  planes 
with  the  free-fleld  stresses  to  be  measured.  He  present  results  for  grout 
fill  in  tuff  that  show  how  this  stress  comparison  is  Influenced  by  the 
following: 

•  Mismatch  of  borehole  fill  and  aedlua  properties 
e  Bonding  of  the  borehole  fill  to  the  aedlua 

•  Separation  of  the  borehole  fill  froa  the  aedlua 
e  Formation  and  distribution  of  residual  stresses. 

Inclusion  effects  caused  by  the  gage  package  are  addressed  in  Section  4. 

He  have  constructed  simple,  practical  exaaples  to  Illustrate  clearly 
the  above  inclusion  effects  on  stress  asasureaent.  To  simplify  the 
example,  we  aade  several  reasonable  assumptions  or  introduced  limita¬ 
tions  on  the  applicability.  Our  exaaples  are  quasi-static  and  apply  to 
common  cases  where  the  rise  time  of  the  spherical  stress  wave  is  greater 
than  the  duration  of  several  wave  transits  across  the  borehole  diameter. 
The  spherical  wave  engulfing  the  borehole  is  assumed  to  generate  a 
stress  field  that  can  be  adequately  represented  by  a  biaxial  stress- 
plane  strain  state,  as  shown  in  Figure  2.  The  coapreeslve  stresses  py 
and  px  represent  the  radial  and  circumferential  coaponents  of  the 
spherical  wave.  Another  simplification  was  to  adopt  proportional 
loading;  that  is,  to  keep  the  stress  ratio  Px/py  constant.  Once  the 
dependency  on  the  proportional  loading  ratio  is  understood,  the  borehole 
stress  histories  can  be  interpreted  froa  more  realistic  loading  paths 
(spherical  wave  coaponents).  The  grout  and  tuff  properties  have  been 
approximated  in  our  initial  calculations  by  elastic,  perfectly  plastic 
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behavior  with  yielding  and  plastic  deformation  according  to  the  Mises 
yield  condition  and  the  associated  flow  rule.  Mohr-Coulomb  yielding  and 
Irreversible  compressibility  will  be  included  in  a  later  phase  of  the 
program.  For  the  grout  chosen,  designated  2C4,  and  for  several  NTS 
tuffs,  the  friction  angles  are  not  large,  so  the  Mlses  criterion  is  a 
reasonable  approximation.  It  is  probably  more  important  to  Include 
plastic  compressibility  in  the  future. 

The  material  properties  used  in  the  calculations  are  listed  in 
Table  1. 


Table  1.  Material  Properties 


Material 


Tuff  medium 


Young’s  Poisson’s 
Modulus,  Patio, 

E  (psl)  v 


1.66  x  10w 


Grout  (2C4)  fill  2.64  x  106  0.28 


Shear 
Modulus , 
u  (psi) 


0.63  x  10* 


1.03  x  10® 


Crush 
Strength, 
o„  (psi) 


Figure  3  shows  the  two  gage  planes  for  measuring  the  free-field 

stress  components  px  and  Py.  Because  of  symmetry,  only  the  first 

quadrant  needs  to  be  modeled.  Stresses  o  and  o  are  calculated  in  the 
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inclusion  along  the  radii  on  the  y  and  x  axes,  respectively.  The 
inclusion  effect  of  the  gage  itself  is  not  considered  in  these  borehole 
calculations,  and  the  stresses  ax  and  oy  along  the  y  and  x  axes  are 
taken  to  be  the  stresses  acting  on  the  two  gages  for  comparison  with  the 
uniform  free-field  stresses. 

The  analyses  were  performed  with  the  finite  element  code  NIKE. 

The  finite  element  grid  was  extended  to  10  times  the  inclusion  radius, 
as  shown  in  Figure  4.  Four  node  plane-strain  quadrilateral  elements 
were  used. 
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Three  inclusion-medium  interface  conditions  are  Investigated. 

These  Include  perfectly  bonded,  free-sliding,  and  frictional  sliding 
conditions.  The  perfectly  bonded  and  free-sliding  interfaces  provide 
bounds  on  the  more  realistic  frictional  sliding  interface.  The  bonded 
Interface  transfers  compressive,  tensile,  and  shearing  loads.  The  free- 
sliding  Interface  transfers  compressive  loads  only.  The  frictional 
sliding  interface  transfers  compressive  loads  and  limits  the  amount  of 
shearing  load  that  can  be  transferred. 

The  perfectly  bonded  interface  requires  no  special  modeling 
technique  because  the  usual  finite  element  node-element  connectivity 
describes  this  condition. 

The  free-  and  frictional  sliding  Interface  capabilities  are  based 
on  a  finite  element  penalty  function  formulation.  The  Interface 
geometry  is  defined  by  a  list  of  nodal  point  numbers  lying  along  both 
sides  of  the  Interface.  One  side  of  the  contact  surface  is  referred  to 
as  the  master  surface,  and  the  other  side  is  referred  to  as  the  slave 
surface.  In  a  region  of  contact  and  sliding  friction,  Interface  springs 
are  used  to  preclude  penetration  and  sliding.  Any  node  that  penetrates 
through  or  slides  away  from  its  respective  contact  surface  causes  a 
linear  Interface  spring  to  be  inserted  into  the  stiffness  matrix;  this 
spring  couples  the  penetrating/ sliding  node  to  two  adjacent  nodes  on  the 
contact  surface.  An  adjustable  penalty  parameter  is  input  to  scale  the 
stiffness  of  the  Interface  springs,  and  a  coefficient  of  friction  is 
input  to  prescribe  the  frictional  behavior  of  the  Interface. 

The  numerical  results  treat  the  effects  on  the  stress  distributions 
acting  normal  to  the  gage  planes  of  the: 

•  Interface  conditions  (tuff  inclusion) 

e  Material  properties  (grout  inclusion) 

•  Combined  interface  conditions  and  material  properties 

•  Unloading  (residual  stress) 

e  Interface  separation. 
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In  all  but  one  case,  the  media  were  taken  to  be  an  elastic- 
perfectly  plastic  Von  Mlaea  nodal  of  tuff*  The  calculations  were  nade 
for  two  or  three  free-fleld  stress  levels,  the  lowest  value  In  each  case 
providing  an  entirely  elastic  response. 

2.2  INTERFACE  CONDITIONS 

The  ideal  borehole  inclusion  would  be  perfectly  bonded  to  and  natch 
the  material  properties  of  the  surrounding  media.  If  this  state  could 
be  achieved,  the  stress  distribution  in  the  borehole  would  be  identical 
to  the  free-fleld  stress. 

Figure  S  shows  the  stress  distributions  in  an  inclusion  that 
matches  the  material  properties  of  the  media  exactly,  but  allows  the 
inclusion  to  slip  freely  along  the  Interface.  The  stress  distributions 
along  the  x  and  y  axes  are  shown  for  a  free-fleld  compressive  stress  of 
Py  ■  6  and  8  ksl,  with  a  proportional  loading  ratio  of  Py  “  4PX.  A  gage 
along  the  x  or  y  axis  would  measure  the  free-fleld  stress  exactly  if  the 
ratios  of  ax/V%  and  oy/Py  are  unity.  The  effect  of  a  slip  Interface 
that  does  not  transfer  shear  stress  is  to  produce  a  nonuniform  stress 
distribution  along  both  axes.  At  the  lower  loading  level,  Py  ■  6  1, 

the  tuff  Inclusion  is  elastic  (solid  lines),  and  at  the  center  of  the 
inclusion,  the  x  and  y  stress  components  are  0.27  and  1.17  of  their 
free-fleld  counterparts,  respectively.  Increasing  the  free-fleld 
loading  to  8  ksl  causes  yielding  in  the  Inclusion  and  flattens  the 
stress  distribution  at  the  center,  where  the  stress  ratios  are  closer  to 
unity  (dashed  lines).  Approximate  measurement  of  the  large  free-fleld 
stress,  Py,  is  possible  if  the  gage  does  not  occupy  more  than  50Z  of  the 
inclusion  diameter.  Measurement  of  the  smaller  free-field  component, 

PR,  can  be  regarded  only  as  an  estimate,  the  accuracy  of  which  depends 
on  the  free-field  stress  level. 

Figure  6  shows  the  stress  distributions  along  the  x  and  y  axes  of  a 
tuff  inclusion  with  a  frictional  Interface  with  a  coefficient  of  p  *  0.6 
for  free-field  compressive  stresses  of  Py  “  4,  6,  and  8  ksl,  and  a  load¬ 
ing  ratio  of  Pv  -  4P_.  The  effect  of  including  a  frictional  Interface 


Figure  5.  Stress  along  gage  planes  for  tuff  inclusion 
with  slip  interface. 


is  to  produce  More  uniform  stress  distributions  than  those  obtained  with 
a  slip  interface.  The  stress  ratios  (meaaured/free-field)  of  the  x  and 
y  components,  respectively,  are  0.66  and  0.98  for  elastic  response  and 
0.78  and  0.95  for  plastic  response  at  Py  •  8  ksl.  The  accuracy  in 
measurement  of  both  free-field  stresses,  Px  and  Py,  Increases  in  going 
from  a  slip  to  a  frictional  interface.  The  accuracy  of  measuring  the 
larger  free-field  stress,  Py,  would  rapidly  approach  that  obtainable  for 
a  bonded  interface  as  the  coefficient  of  friction  is  increased. 

The  coefficient  of  friction,  p  *  0.6,  used  in  all  the  frictional 
interface  calculations  is  an  estimate  of  the  sliding  friction  coefficient 
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for  tuff  on  tuff,  based  on  a  survey  paper  by  Byerlee.  Methods  for 
measuring  the  friction  coefficients  of  tuff  on  tuff  and  tuff  on  grout 
for  both  sliding  friction  and  bonding  strength  of  a  poured  Interface 
will  be  described  in  a  future  report. 

2.3  MATERIAL  PROPERTIES 

In  this  section,  we  investigate  the  effect  of  having  a  borehole 
with  material  properties  that  do  not  match  the  surrounding  media.  The 
mismatch  of  a  strong  tuff  and  a  weak  grout  (2C4)  is  intentionally  chosen 
to  magnify  the  effects  of  poorly  matched  properties. 

Figure  7  shows  uniform  stress  distributions  along  the  x  and  y  axes 
of  a  bonded  grout  inclusion  for  free-field  compressive  stresses  of  Py  ■ 

4,  6,  and  8  ksl,  with  a  proportional  loading  ratio  of  P  «  4P  .  The 
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stress  ratio  for  the  larger  free-field  component,  Py»  starts  at  1.1  for 
elastic  response  and  decreases  to  0.93  for  plastic  response  of  the 
inclusion  at  P  •  8  ksl.  Conversely,  the  smaller  free-field  stress,  P  , 
starts  at  0.93  for  elastic  response  and  Increases  to  1.49  for  plastic 
response  at  Py  ■  8  ksl.  Because  all  the  points  in  the  inclusion  are  on 
the  yield  surface,  this  trend  in  stress  ratios  would  continue  with 
increasing  load  until  a  substantial  portion  of  the  medium  near  the 
borehole  also  reached  the  yield  surface.  When  a  significant  portion  of 
the  medium  yields,  there  would  be  a  redistribution  of  load  near  the 
inclusion,  which  could  change  the  stresses  in  the  borehole.  Even  with 


th«  mismatch  of  aaterial  properties  la  this  exsap le,  s  reasonably 
accurate  aaasureaent  of  the  larger  free-field  stress,  Py,  is  obtained. 
The  seas  is  true  for  aeesureaent  of  the  saaller  coaponent  under  elastic 
or  low  level  plastic  loading.  As  the  loading  Increases  to  produce  aore 
plastic  deforest ions,  asasureaents  of  Px  becoae  aore  Inaccurate. 

2.4  INTERFACE  CONDITIONS  AND  MATERIAL  PROPERTIES 

Figure  8  shows  the  stress  distributions  along  x  and  y  axes  of  a 
grout  Inclusion  with  a  slip  Interface  for  free-field  coapressive 
stresses  of  Py  -  1.5,  4,  and  6  kal,  with  Py  -  4PZ.  This  exaaple 
represents  the  worst  case  considered  because  it  coablnes  the  effects  of 
aisaatched  aaterial  properties  and  the  lack  of  Interface  bonding.  The 
effect  of  aaterial  aisaatch  is  only  significant  in  aeasureaents  of  the 
saaller  free-field  streas,  ?x,  as  seen  in  Figure  7.  Cowparing  Figures  7 
(bonded  grout  Inclusion)  and  Figure  8  (slip  grout  inclusion),  we  con¬ 
cluded  that  the  effect  of  bonding  la  to  produce  unifora  stress  distribu¬ 
tions  and  to  Increase  the  accuracy  of  strass  aeasureaents  at  the  lower 
free-field  load  levels.  These  conclusions  are  also  illustrated  in 
Figure  9,  which  shows  the  stress  distributions  along  the  x  and  y  axes  of 
a  grout  inclusion  with  a  frictional  interface.  The  results  for  the 
frictional  Interface  are  bounded  by  the  bonded  and  slip  Interface 
results  as  expected. 

Based  on  the  nuaerlcal  exaaples  for  loading  presented  in  Sections 
2.2,  2.3,  and  2.4,  we  conclude  the  following: 

•  Borehole  aaterial  should  be  well  bonded  to  the  aedlua  to 
ensure  a  uni fora  stress  distribution. 

•  Borehole  asterlsl  properties  should  aatch  the  aedla 
properties  for  iaproved  sccuracy. 

•  Multiple  stress  coaponent  aeesureaent  is  difficult  to 
obtain  directly. 

a  Ratio  of  borehole  stress  coaponent  to  corresponding  free 
field  is  not  constant  if  slip  occurs. 

•  Uniform  (±102)  stress  regions  exist  for  <502  of  the  plsne 
of  s  gage. 
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Figure  8.  Stmt  along  gage  planet  for  grout  inclution 
with  tlip  interface. 
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Figure  9. 


Stress  along  gage  planes  for  grout  inclusion 
with  frictional  interface  p  -  0.6. 


2.5  UNLOADING  AND  RESIDUAL  STRESS 


When  the  material  In  the  borehole  Inclusion  flows  plastically 
before  and  during  unloading,  the  mismatching  of  stress  and  displacement 
across  the  borehole  Interface  and  the  mismatching  of  borehole  material 
properties  will  result  In  a  residual  stress  field.  This  residual  stress 
field  can  prevent  gage  signals  from  returning  to  their  original  baseline 
and  cauae  inaccurate  impulse  measurements.  In  this  section,  we 
Investigate  the  magnitude  of  this  problem. 

The  residual  stress  field  is  dependent  on  the  entire  history  of 
loading  and  unloading.  Including  proportional/nonproportional  loading, 
maximum  stress  attained  during  loading,  material  characterization 
(constitutive  model),  and  proportional /nonproportional  unloading.  To 
Illustrate  these  dependencies,  we  considered  a  bonded  grout  inclusion 
surrounded  by  tuff  that  Is  loaded  and  unloaded  proportionally. 

Figure  10  shows  the  stress  space  trajectories  of  any  material  point 
In  a  bonded  grout  inclusion  using  Tresca  and  Von  Mises  yield  criteria 
with  two  proportional  loading  and  unloading  ratios.  Figure  10(a)  shows 
the  stress  trajectories  for  the  inclusion  with  a  proportional  loading 
and  unloading  ratio  of  Py/Px  ■  •  (l.e.,  Px  “  0)  for  Tresca  (dashed  line) 
and  Von  Mises  (solid  line)  yield  criteria.  Figure  10(b)  uses  the  same 
notation  for  the  proportional  loading  ratio  Py/Px  "A. 

The  Von  Mises  yield  stress  trajectory  passes  through  the  Tresca 
yield  surface  until  It  Intersects  the  Von  Mises  yield  surface.  (Note: 

In  the  x-plane,  the  Von  Mises  yield  surface  is  a  circle  that  circum¬ 
scribes  the  Tresca  regular  hexagon).  Yielding  continues  along  each  of 
the  yield  surfaces  until  the  larger  free-field  stress,  Py,  reaches 
8  kel.  The  tuff  media  are  then  unloaded  using  the  same  ratio  of  P„/P_ 

J  * 

as  In  the  loading  phase.  The  unloading  curves  are  parallel  to  the 
original  elastic  loading  curves  unless  reverse  yielding  occurs  on 
unloading,  as  it  does  along  the  Tresca  trajectory  in  Figure  10(a). 

The  beginning  and  ending  of  the  unloading  trajectories  are 
indicated  by  X's  for  the  Von  Mises  yield  criteria  and  squares  for  the 


Tresca  yield  criteria  in  Figure  10.  The  final  stress  states  far  the  two 
Von  Nises  trajectories  are  1275  psi  apart;  the  Tresca  points  are  582  psi 
apart.  Coaparlson  of  the  final  stress  states  shows  that  in  Figure  10(a) 
for  the  loading/unloading  ratio  of  Py/Px  "  •»  the  Von  Mises  and  Tresca 
states  are  308  psi  apart,  and  in  Figure  10(b)  for  Py/Px  “  4,  the  two 
points  are  447  psi  apart.  Froa  these  results,  we  conclude  that  residual 
stress  distributions  are  more  dependent  on  loading  and  unloading  ratios 
than  constitutive  relation.  The  choice  of  constitutive  relation  is 
probably  more  significant  for  rocks  with  higher  friction  angles  where 
Mohr-Coulomb  type  yield  criteria  are  appropriate  and  the  mean  pressure 
affects  the  residual  stress  field. 

Effect  of  Interface  Conditions  on  Residual  Stress 

Figure  11  shows  the  residual  stress  distributions  in  a  tuff 
inclusion  that  was  loaded  to  a  f ree-f ield  stress  of  Py  ■  8  Itsl  with  a 
loading  ratio  of  Py/Px  “  4.  The  residual  stress  distributions  are 
plotted  for  a  slip  interface  (solid  lines)  and  a  frictional  interface 
(dashed  lines).  The  bonded  interface  does  not  result  in  a  residual 
stress  distribution.  The  resulting  residual  stress  distributions  have 
the  same  dependency  on  the  interface  condition  that  was  found  in  the 
loading  examples.  Frictional  interfaces  produce  low  level  uniform 
residual  stress  distributions,  whereas  slip  interfaces  produce  large 
residual  stresses  that  depend  on  location  in  the  inclusion.  The 
difference  in  residual  stress  distributions  for  the  frictional  and  slip 
interfaces  is  that  for  bonding  or  high  friction,  all  material  points  in 
the  inclusion  are  at  the  same  stress  state.  In  the  frictional  interface 
inclusion,  all  points  also  unload  uniformly  and  distribute  the  residual 
stress  uniformly.  The  slip  interface  causes  a  distinct  elastic-plastic 
zone  to  develop  within  the  inclusion  during  loading,  and  hence  causes 
differing  amounts  of  residual  stress  during  unloading. 

Figure  12  shows  the  residual  stress  distributions  for  a  grout 
inclusion  that  has  been  loaded  to  a  maximum  free-fleld  stress  of  Py  - 
6  ksl  with  Py  ■  4PX«  The  residual  stress  distributions  for  three 
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Figure  12.  Residual  stress  along  gage  planes  with  grout  inclusion. 

interface  conditions  are  plotted:  bonded  interface  (solid  lines),  slip 
interface  (short  dashes),  and  frictional  interface  (long  dashes).  The 
bonded  interface  results  in  a  large  uniform  residual  stress  distribution 
The  slip  interface  results  in  a  residual  stress  field  with  the  same 
distribution  as  that  calculated  for  the  tuff  inclusion.  The  frictional 
Interface  produces  the  lowest  level  of  residual  stress  and  has  a 
distribution  similar  to  the  slip  interface. 

The  results  presented  in  this  section  confirm  some  of  the 
conclusions  presented  for  the  loading  examples: 

e  Borehole  material  should  be  bonded  to  the  media  to  ensure 
uniform  stress  distributions. 

•  Borehole  material  properties  should  match  the  media 
properties  for  minimal  residual  stress. 

e  Multiple  stress  component  measurement  is  difficult  to 
achieve  directly. 

2.6  INTERFACE  SEPARATION 
Interface  Separation  Formulation 

Clearly,  stress  measurements  will  be  affected  if  the  inclusion 
separates  from  the  media.  When  elastic  conditions  prevail,  the  ratios 
of  the  free- fie Id  stress  components  that  cause  radial  tension  across  the 
incluslon/medlum  interface  can  be  determined  explicitly.  In  the  more 
general  case  of  elastic-plastic  loading  and  unloading,  numerical 
algorithms  will  be  needed  to  determine  the  stress  distributions.  The 
numerical  separation  algorithm  was  verified  by  simulating  a  sample 
problem  whose  analytical  solution  is  known. 

The  example  considers  the  plane-strain  elastic  contact  between  a 
circular  disk  and  a  circular  hole  in  an  infinite  plate.  The  disk  is 
smooth,  and  the  diameter  of  the  disk  and  the  hole  are  identical  in  the 
unstressed  state.  A  uniaxial  stress  is  applied  to  the  plate  at 
infinity,  resulting  in  a  separation  along  a  portion  of  the  boundary 
between  the  plate  and  disk  as  indicated  in  Figure  13(a).  The  analytical 
solution  for  the  arc  of  contact  and  the  stress  distribution  along  this 
arc  are  given  in  Reference  4. 

37 


Figure  13(b)  shove  the  radial  stress  distribution  along  the  contact 
arc,  an  angle  of  19.62*,  for  the  analytical  solution  (solid  line)  and 
the  calculated  values  in  the  plate  (+'s)  and  in  the  inclusion  (X's). 

The  numerical  results  are  in  good  agreenent  with  the  analysis  and 
probably  could  be  improved  by  refining  the  computational  mesh. 

Figure  13(c)  shorn  analytical  solutions  of  the  hoop  stress  around 
the  hole  in  an  infinite  plate  for  the  contact  of  a  smooth  disk  with  a 
circular  hole  in  an  infinite  plate  (solid  line);  for  comparison,  an 
Infinite  plate  with  a  circular  hole — that  is,  no  disk  (dashed  line)  is 
shown.  The  computed  hoop  stress  for  the  contacting  disk  (X's)  shown  in 
Figure  13(c)  is  in  very  good  agreement  with  the  analysis  of  the 
contacting  disk. 

Stress  Distributions 

If  we  let  the  uniaxial  free-field  stress  be  compressive  so  that 
■  -p^  (p^  >  0,  «*  0),  as  shown  in  Figure  14  (lower  left),  the 

contact  regions  are  centered  on  the  y  axis.  The  compression  py  was 
chosen  at  a  value  of  1000  stress  units.  As  long  as  the  material  remains 
elastic,  the  stress  distributions  remain  the  same  relative  to  this  free- 
fleld  compressive  stress.  Figure  14  shows  the  stress  distributions  oy 
and  ox  along  the  planes  y  ■  0  and  x  •  0,  that  is,  along  the  gage 
planes.  The  stress  distributions  are  nonuniform,  being  qualitatively 
parabolic  in  shape.  The  main  stress  component  of  Interest,  Oy,  is  Just 
over  20Z  high  near  the  center  of  the  borehole  fill,  but  stress  measure¬ 
ment  is  still  useful  if  the  gage  is  confined  to  the  central  half  of  the 
fill  diameter. 


Hi 


SECTION  3 


STATIC  ANALYSIS  OF  INCLUSIONS 

3.1  INTRODUCTION 

/ 

We  have  obtained  formulas  for  the  stresses  in  elliptic  and  circular 
Inclusions  under  plane  strain  conditions.  The  results  apply  when  the 
Inclusion  and  medium  remain  elastic,  and  we  obtained  them  by  performing 
analyses  or  adapting  known  analyses  to  solve  our  particular  problems. 

The  purposes  of  the  analytic  treatment  are  to  validate  the  finite 
element  code  and  to  facilitate  parametric  investigations  In  the  elastic 
regime.  The  elliptic  inclusion  results  assist  the  study  of  the  Instal¬ 
lation  Inclusion  formed  when  the  gages  are  inserted  In  slots  from  a 
trench  wall.  These  results  are  also  useful  in  developing  stress  gages 
of  ytterbium  foil  In  shear  compliant  material  for  laboratory  applica¬ 
tions,  such  as  in  DNA  containment  research.**  The  circular  inclusion 
results  apply  to  the  Inclusion  formed  by  borehole  fill. 

The  specific  Inclusion  problems  addressed  are: 

e  Bonded  general  ellipse 

•  Bonded  slender  ellipse 
e  Bonded  circle 

•  Unbonded  circle  (no  Interface  shear) 
e  Separation  conditions 

e  Residual  stress  in  bonded  circle  (elastic-plastic  inclusion) 
e  Rigid  line  inclusion. 

The  complex  variable  method  of  Muskhellshvlli^  Is  used  to  obtain  the 
solutions,  except  when  the  circular  inclusion  becomes  plastic  during  the 
residual  stress  investigation. 


3.2  BONDED  ELLIPTIC  INCLUSIONS 


Figure  15  shove  an  elliptic  Inclusion  under  plane  strain  conditions 
In  a  aedlua  with  far-fleld  principal  stresses  ax  m  ~VX»  ay  m  ~P^ •  For  a 
perfectly  bonded  interface,  the  stress  distribution  in  the  inclusion  is 
uniform  and  equal  to 

-  “i-h 


where 


a  k2[{a  k2  -  (a  -  1)<1  +  a2)}^  -  (a  -  l)a  Gj 


kj[a  k j  -  (a  -  1)(1  -  ■  )]  +  2(«  -  1)(1  -  ■*)[«  kj  -  <«  -  1)] 
-  [kt  +  2(a  -  1)1A1  +  «  k2  Ct 


"  T  <py  +  G1  “  "  T  (Py  “  Vx> 


The  subscript  1  refers  to  properties  within  the  Inclusion,  andthe 
subscript  2  refers  to  properties  within  the  aedlua.  Kj  and  K2  are  the 
bulk  aodull,  and  p2  are  the  shear  aoduli,  and  and  v2  are  Poisson's 
ratios. 

We  use  the  shore  result  to  deteraine  the  stresses  in  a  bonded 
slender  ellipse  and  in  a  bonded  circle  for  the  study  of  a  gage  package 
and  a  borehole  fill  Inclusion,  respectively. 


3.3  BONDED  SLENDER  ELLIPTIC  INCLUSION 


The  slender  elliptic  inclusion  is  of  interest  because  it  represents 
an  overall  elastic  behavior  of  a  flat pack  and  the  elastic  behavior  of  a 
laboratory  stress  gage  with  an  ytterbium  foil  embedded  in  fluid  or  a 
shear  compliant  notarial  such  as  silicone. 


Let  the  ratio  of  the  ainor  and  aajor  axes  be  b/e  "  e,  where  e  is  a 
small  number  defining  the  slenderness  ratio.  Substitute  a  •  (1  -  e)/ 

(1  +  e)  in  Equation  (1)  and  neglect  terns  in  e  higher  than  the  first 
order.  Then  we  obtain  for  the  stress  oy  the  formula 


The  stress  within  the  inclusion  is  dependent  upon  the  aspect  ratio  of 
the  ellipse,  e  “  b/a;  the  aaterlal  properties,  v^,  v^,  and  a  “ 
and  the  free-fleld  stress  ratio  px/py.  When  the  material  in  the 
inclusion  is  the  same  as  the  medium  (a  ■  1,  kj  ■  k^) ,  Equation  (2)  gives 
a  »  -p  ,  as  it  should. 

y  y 

As  the  aspect  ratio  of  the  ellipse  becomes  small  (c  +  0),  Equation 
(2)  gives  oy  ♦  ~Py .  If  the  ellipse  can  be  made  very  slender,  it  does 
not  aatter  whether  or  not  the  material  elastic  properties  in  the  inclu¬ 


sion  natch  those  in  the  medium, 
varies  with  the  property  ratio  a 
ellipse  of  aspect  ratio  e  “  0.1. 


figure  16  shows  how  the  ratio  -a  / p 

y  y 

■  and  loading  ratio  Px/py  for  an 

As  in  the  case  of  grout  and  tuff,  we 


assume  that  the  values  of  the  Poisson's  ratio,  v,  are  comparable,  so  for 
simplicity  we  have  chosen  v  “  0.3,  the  average  of  the  values  in  Table  1. 
Over  the  range  1/2  <  ^/^  <  2,  whenever  |px/py|  <1)W  aee  that  the 

Inclusion  stress  coaponent  -cy  is  within  5Z  of  the  free-field  stress 
component  py.  When  the  inclusion  aaterlal  is  weaker  than  the  mediun 
(1/2  <  P1/112  ^  the  inclusion  stress  is  lower  than  the  free-field 
stress.  When  the  inclusion  aaterlal  is  stronger  (1  <  / M-2 ) »  the 

Inclusion  stress  is  higher  than  the  free-field  stress. 


In 


summary,  if  we  wish  to  measure  py,  the  slender  elliptic 
Inclusion  allows  us  to  do  this  with  acceptable  accuracy  unless  the 
lateral  component  px  becoaes  too  large,  say  |px/py|  >  2. 


In  one  design  of  stress  gage,  the  ytterbium  foil  is  surrounded  by  a 
material  that  may  be  considered  a  fluid,  that  is,  0  and  -  1/2 
a  *  0,  »  2).  The  pressure  in  the  fluid  Inclusion  is  determined  by 

Equation  (2)  as 


-L  i  +  1  ~  2v2  f.  3  K2 

Py  "  6  1  •  v2  l  2(1  +  v2}  Kl. 


(3) 


Equation  (3),  to  the  first  order  in  e,  shows  that  the  inclusion  pressure 
is  independent  of  the  free-field  stress  component,  px«  Figure  17  shows 
the  variation  of  the  inclusion  pressure  with  the  medium  Poisson's  ratio, 
V£,  and  the  ratio  of  the  bulk  moduli,  K^/Kj*  The  curves  show  that  the 
inclusion  pressure  is  within  5Z  of  the  free-field  stress,  py,  when  the 
fluid  bulk  moduli  is  greater  than  the  medium  bulk  modulus  by  up  to  a 
factor  of  three.  They  also  show  that  the  fluid  modulus  should  not  be 
leas  than  half  the  medium  bulk  modulus,  based  on  about  a  90X  accuracy 
when  V2  *  0.3. 


3.4  BONDED  CIRCULAR  INCLUSION 

The  stress  oy  in  a  bonded  circular  inclusion  representing  the 
borehole  fill  is  obtained  by  specializing  Equation  (1)  to  the  case 
m  ■  0,  giving  a  ■  b.  This  process  leads  to 


and  where  again  a  ■  / 1*2 *  ^1  “  ”  v),  k2  ■  4(1  -  V2>*  The 

degree  to  which  the  stress,  y,  measured  within  the  inclusion  matches 
the  free-f ield  stress  component,  Py,  depends  on  the  ratio  a  of  the  shear 
moduli,  Poisson’s  ratio  and  V2«  end  the  ratio  Pz/Py  °f  the  free-field 
8 tress  components*  For  the  special  case  in  which  the  inclusion  and 
medium  are  the  same  material  (a  ■  1,  lj_  ■  k2)»  Equation  (4)  gives 
oy  ■  -py,  as  it  should. 

Figure  18  shows  how  the  ratio  -<Ty/py  varies  with  the  property  and 
loading  ratios  p^p^  and  Px/Py*  As  In  the  case  of  grout  and  tuff,  we 
assume  that  the  values  of  the  Poisson's  ratio,  v,  are  comparable,  so  for 
simplicity  we  have  chosen  v  ■  0.3,  the  average  of  the  values  in  Table  1. 
The  curves  pass  through  the  point  (1,1)  representing  the  case  of 
Identical  inclusion  and  medium  materials  (or  no  Inclusion).  When  the 
Inclusion  is  stiff er  than  the  medium  (Pj/Pj  *  1)»  the  stress  measured  in 
the  Inclusion  is  high  (-  cry/py  >  1).  For  p^  •  2p^  and  px  “  0,  the 
measurement  is  about  20Z  high;  the  accuracy  Improves  as  px  Increases. 
When  the  inclusion  is  more  compliant  than  the  medium  <  1),  the 

stress  measured  in  the  inclusion  is  low  (-  Oy/py  ^  1) •  For  p^  ■  P2/2 
and  px  "  0,  the  measurement  is  about  25Z  low;  again,  the  accuracy 
Improves  as  px  Increases. 

Figure  18  illustrates  that  when  elastic  behavior  dominates,  it  is 
good  practice  to  match  the  medium  modulus  with  an  Inclusion  modulus  to 
within  about  20Z  if  lower  than  the  medium  modulus  and  about  40Z  if 
higher. 

3.5  UNBONDED  CIRCULAR  INCLUSION 

If  the  gage  installation  procedure  results  in  a  bond  between  the 
borehole  fill  Inclusion  and  the  medium  that  is  weak  in  shear  strength, 
the  loading  may  cause  slip  to  occur  along  the  interface.  The  slip 
radically  alters  the  stress  distribution  in  the  Inclusion.  In  the 
extreme  case  of  an  interface  with  no  shear  strength,  this  stress 
distribution  is 


Figure  18.  Stress  in  a  bonded  circular  inclusion  (»>.  =  v7  «  0.3). 
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where 


P  *  T  1<«  -  1)  and  6 


kL  +  3ak2  "  2(o  -  1) 


a  ■  ki  "  *<1  ~  vt)  k2  -  4(1  -  v2) 

Equation  (5)  shows  that  the  Inclusion  stress  depends  on  the  material 

properties  through  Poisson's  ratios  and  V2  and  the  ratio  a  of  the 
shear  moduli,  and  on  the  loading  ratio  Px/py»  as  in  the  case  of  the 
bonded  Interface.  However,  in  stress  distribution  of  each  component  is 
no  longer  uniform,  but  is  parabolic  in  the  Cartesian  coordinate  system. 

A  stress  gage  located  on  the  x  axis  for  measuring  the  stress  oy 

would  be  subjected  to  a  parabolic  stress  distribution  with  a  maximum 

compressive  stress  at  the  center  (x  »  0,  y  -  0)  when  py  >  p%.  (Note 
that  the  parameters  0  and  6  are  positive.)  For  the  same  case,  a  stress 
gage  located  on  the  y  axis  for  measuring  the  stress  ox  would  be 
subjected  to  a  parabolic  stress  distribution  with  a  minimum  compressive 
stress  at  the  center. 

Figure  19  shows  the  variation  of  -oy/py  with  the  shear  modulus 
ratio,  a,  and  the  loading  ratio  px/py  when  “  v2  *  0*3 »  a  value 
appropriate  for  NTS  tuff.  The  stress  values  are  for  the  center  of  the 
inclusion;  we  set  x  ■  0  in  Equation  (5).  As  stated  above,  if  p_/p_  <  1, 
Figure  19  shows  the  maximum  compressive  stress,  and  if  Px/py  >  1,  Figure 
19  shows  the  minimum  compressive  stress.  If  p  *  p  ,  the  stress 
distribution  is,  ot  course,  uniform.  If  the  gage  occupies  a  small 
enough  portion  of  the  diameter  on  the  x  axis,  the  influence  of  the 
parabolic  distribution  will  be  reduced  and,  as  an  approximation,  the 
gage  stress  will  be  determined  by  Figure  19.  We  see  immediately  that 
stress  measurement  is  reasonable  only  when  py  >  px,  although  estimates 
are  possible  for  values  of  px  in  the  range  pv  <  p,  <  3p  . 
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unbonded  circular  inclusion  (»>.  =  v0  =  0.3;  x  =  0). 
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As  a  final  comment  on  Figure  19 ,  when  values  of  the  loading  ratio 
px/py  are  less  than  about  0.5  and  greater  than  about  5  (depending  on  the 
shear  modulus  ratio),  the  inclusion  and  medium  can  separate  if  the 
Interface  bond  has  negligible  tensile  strength.  Separation  will  affect 
the  stress  distribution.  The  upper  curve  indicates  that  separation  will 
occur  at  0  ■  0  for  low  values  of  the  lateral  loading  px»  Separation 
will  occur  at  0  ■  90°  for  high  values  of  the  lateral  stress  when  the 
px/py  curve  intersects  the  abscissa,  the  a  axis. 

Figure  20  shows  the  variations  of  -oy/py  with  the  loading  ratio, 
px/py,  at  several  points  along  the  positive  x  axis  in  the  inclusion  when 
the  inclusion  and  medium  are  2C4  grout  and  tuff  with  the  elastic  proper¬ 
ties  listed  in  Table  1  "  1.63).  At  each  radius,  given  by  x/a, 

the  curves  are  straight  lines  in  accordance  with  Equation  (5).  The 
farther  the  values  of  the  loading  ratio  are  from  unity,  the  more  pro¬ 
nounced  is  the  parabolic  stress  distribution  along  the  x  axis,  and  the 
parabola  inverts  as  we  pass  through  px/py  ■  1» 

Figure  21  shows  the  variation  of  ~cry/py  with  the  loading  ratio, 
px/Py,  at  the  center  and  mldradlus  of  an  inclusion  of  2C4  grout  in  a 
tuff  medium,  again  frith  the  elastic  properties  listed  in  Table  1.  The 
third  line  (dashed)  is  bonded  interface  according  to  Equation  (4);  for 
the  bonded  case,  the  stress  oy  applies  throughout  fhe  entire  inclusion. 
The  lines  indicate  the  importance  of  bonding  when  measuring  stresses  in 
the  elastic  regime,  as  indicated  by  the  departures  from  the  desired 


value  of  unity  for  -cy/py. 


3.6  SEPARATION  AT  THE  INTERFACE 


By  examining  Equation  (5)  for  the  unbonded  Interface,  we  have  seen 
that  separation  can  occur  at  the  interface.  In  this  section,  we  present 
the  separation  condltons  for  the  bonded  and  unbonded  Interfaces  with  no 
tensile  strength. 


The  radial  component  of  stress  is 


Jr  "  “  7  (px  +  py>  0  "  7  (px  "  py5  6  co*  20 


I 


Stress  in  unbonded  2C4  grout  circular  inclusion  in  tuff  (properties  in  Table  t; 

V\l»2  *  163>- 


+  2(0  -  lj 


-  (a  - 


(bonded  or  unbonded) 


(bonded) 


-  21  a  - 


(unbonded) 


where  a  -  \  “  v^,  and  k2  -  4(1  -  v2).  Separation  occurs 

et  0  *  0*,  u  when 

h.  •>  «£1  (7a) 


end  at  8  ■  «/2,  3w/2  when 


p  '  6  +  0 


To  Illustrate  the  Influence  of  the  Shear  modulus  ratio  a,  we  have  chosen 
»  v2  •  0.3  and  have  plotted  Equations  (7a)  and  (7b)  as  equalities  in 
Figure  22(a)-(d) .  The  curves  then  give  the  values  of  the  loading  ratio 
py/px  that  cauaa  separation. 

3.7  RESIDUAL  STRESS  IN  A  BONDED  CIRCULAR  INCLUSION 

When  the  stresses  are  high  enough  to  cause  plastic  deformation  in 
the  inclusion,  stresses  reuain  in  and  around  the  inclusion  after 
unloading.  In  the  case  of  the  borehole  fill  inclusion,  a  stress  gage 
will  therefore  record  stress  even  after  the  loading  from  the  stress  wave 
has  passed.  If  the  residual  stress  acting  on  a  flatpack  is  compressive, 
the  signal  will  not  return  to  the  baseline. 

We  have  performed  a  residual  stress  analysis  for  the  case  of  an 
elliptic  inclusion  under  plane  strain  conditions  subject  to  the 


No  Separation 


Bonded  circular  inclusion  separation  conditions  at  0  =  0°  (p.  ■  *»-  *  0.3). 


UVV*. 


Figure  22d.  Bonded  circular  inclusion  separation  conditions  at  6  *  90°  (»».  =  v~  =  0.3). 


following  simplifying  assumptions: 

•  The  Inclusion  Is  elastic-plastic,  whereas  the  medium 
remains  elastic* 

•  The  plastic  behavior  is  governed  by  the  Mohr-Coulomb  yield 
condition  and  Its  associated  flow  rule. 

•  The  loading  Is  proportional  (px  -  Xpy)» 

•  The  inclusion/medium  interface  is  bonded. 

We  illustrate  the  results  by  confining  ourselves  to  the  case  of  a 
circular  borehole  fill  inclusion  with  yielding  of  the  grout  in  accor¬ 
dance  with  the  Tresca  condition  and  the  associated  flow  rule. 

The  elastic  loading  creates  the  Inclusion  stress  field 


or  -  [1  (1  +  X)  0  +  (1  -  X)  6]  p  12 
x  7 

or  -  [-  (1  +  X)  0  -  (1  -  X)  6]  p  12  (8) 

7  7 

where  X  “  px/py  1®  the  far-field  loading  ratio,  and  the  elasticity 
parameters  0  and  6  are  those  of  the  bonded  circular  inclusion  given  in 
Equation  (4).  For  X  <  1,  the  appropriate  Tresca  criterion  is 


(9) 


where  oc  >  0  is  the  grout  crush  strength.  Substituting  the  stresses  of 
Equation  (8)  in  Equation  (9)  gives  the  yield  loading 


Py  -  <JC/(1  -  X)  6  (10) 

We  have  shown  that  further  loading,  which  causes  plastic  flow,  creates 
the  Inclusion  stress  field 

«,--Jo  +  x)  epy  +  £.c 

°y  "  "  7  (1  +  X>  Ppy  "  7  °c  (11) 
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in  which  Py  >  Py.  Let  the  msxlmum  loading  be  Py.  Substituting  of  Py 

for  p  in  Equation  (11)  gives  the  inclusion  stress  field  (a  ,  a  )  at 

*  x  y 

maximum  loading*  Elastic  unloading  by  an  aaount  gives  the  stress 

reduction  (o',  o')  given  directly  by  Equation  (8).  The  stresses  when 
*  _  _ 
the  loading  is  p  -  p'  are  therefore  a  ~  o'  and  a  -  o' .  If  complete 

3^  3^  y  y 

unloading  can  be  performed  elastically,  we  set  -  Py  in  Equation  (8) 
and  subtract  the  result  from  Equation  (11)  to  get  the  residual  stresses 


i 

*  Mp,  2 


p,  <  P,  <  2py 


The  upper  bound  2p  on  the  maximum  loading  expressed  in  Equation  (12)  is 

y 

required  to  aboid  reverse  yielding  according  to  the  criterion 


o-o-o  (1. 

y  x  c 

Substitution  of  Equation  (12)  in  Equation  (13)  shows  that  reverse 
yielding  occurs  when  p^  >  2py. 

If  the  maximum  loading  is  greater  than  twice  the  yield  loading, 
Py  >  2py,  the  stresses  after  elastic  unloading  are 

Ox  -  -  y  (i  +  X)  P  (Py  -  2py)  -  yO{ 

«y  ■  “  J  (I  +  x)  P  (Py  -  2py)  +  y  Oc 

Complete  unloading  then  leads  to  the  simple  result 


3y  "I 


2p  <  p 

y  y 


The  residual  stresses  of  Equations  (12)  and  (14)  would  thus  leave 
two  stress  gages  measuring  px  and  Py  in  compression  and  tension, 
respectively.  The  results  of  this  analysis  were  used  for  comparison 
with  our  code  results  in  Figure  10  (Section  2). 


61 


SECTION  4 


STATIC  CALCULATIONS  OF  STRESS  GAGE  INCLUSIONS 

4.1  INTRODUCTION 

In  this  section,  we  present  the  results  of  a  finite  element 
calculation  for  a  typical  flatpack  stress  gage  embedded  In  2C4  grout. 

The  purpose  of  the  calculation  Is  to  compare  the  normal  stress  experi¬ 
enced  by  the  ytterbium  sensing  foil  to  the  free-field  stress  to  be 
measured.  Following  these  calculations  Is  the  result  of  a  basic  elastic 
analysis  In  which  the  flatpack  Is  Idealized  as  a  rigid  line.  The 
purpose  of  the  analysis  Is  to  examine  the  Influence  of  the  bonding  on 
the  stress  normal  to  the  line. 

The  results  Indicate  excellent  performance  by  the  flatpack  under 
the  static  loading.  If  the  ytterbium  foil  is  confined  to  the  central  60Z 
of  the  flatpack  midplane  width.  The  mldplane  stress  was  not  sensitive 
to  the  bonding  condition  at  the  steel/grout  Interface.  The  idealized 
line  gage  analysis  Indicated  that  the  unbonded  condition  provided  more 
accurate  stress  measurement. 

4.2  FLATPACK  STRESS  GAGE  INCLUSION 

A  cross  section  through  a  typical  flatpack  stress  gage  is  shown  in 
Figure  23.  The  ytterbium  foil  sensing  element  Is  not  shown,  but  it 
occupies  the  mldplane.  We  performed  a  finite  element  calculation  with 
the  flatpack  embedded  in  2C4  grout  and  subjected  to  free-field  stresses 
•  -p^  and  <3y  ”  -p^  of  increasing  levels  while  keeping  the  ratio 
constant  at  Py/px  ■  4.  Plane  strain  conditions  are  imposed.  The  result 
of  Interest  is  a  comparison  of  the  stress  Oy  at  the  mldplane  In  the 
plastic,  which  Is  the  stress  the  ytterbium  foil  measures,  and  that  of 
the  free-field  compressive  stress  py»  Figure  24  shows  the  mldplane 
stress  distributions  for  three  different  free-field  stress  levels — A,  B, 
and  C.  Accuracy  Is  obtained  If  the  ytterbium  foil  does  not  occupy  more 
than  the  central  60Z  of  the  mldplane. 


4.3  LINE  INCLUSION 


As  an  initial  investigation  into  the  effect  on  the  measured  stress 
of  the  bonding  of  a  flatpack  stress  gage  with  the  surrounding  medium,  we 
have  Idealized  the  thin  steel  strip  in  two  different  ways.  In  both 
cases,  we  have  replaced  the  gage  cross  section  by  a  line  of  finite 
length.  In  one  case,  the  line  is  unbonded  to  the  medium  and  supports  no 
shear  traction  along  the  surface.  In  the  other  case,  the  line  is 
lnextensible  and  bonded  to  the  medium. 

Figure  25  shows  the  two  line  problems  that  we  solved  using  the 
method  of  linear  relationship  described  by  Muskhelishvili.®  Only  the 
solutions  for  the  stress,  o^,  normal  to  the  line  will  be  given  here. 

When  the  line  is  unbonded  and  we  essentially  have  a  slot  under 
compression  that  sustains  no  shear  traction,  we  obtain  the  stress  normal 
to  the  line  as 

oy  "  “  7  (PL  +  P2)  +  (Pi  -  P2)  cos  26  (15) 

where  0  is  the  angle  between  the  direction  of  -pj^  in  the  first  quadrant 
and  the  x  axis  (see  Figure  25).  When  0*0,  -  -  p2  "  ~  Py*  E0uation 

(15)  is  exactly  the  y  component  of  stress  in  the  far  field;  hence  if  the 
gage  is  a  perfectly  lubricated  line,  we  obtain  the  correct  reading. 

When  the  line  is  lnextensible  and  bonded  to  the  medium,  we  obtain 
the  stress  normal  to  the  line  as 

ay  “  •  £2Ti  r  (pi  +  p2} "  i  <pi '  p2>  cos  203  {16) 

where  <  •  3  -  4v,  v  being  Poisson's  ratio  for  the  medium.  When  0-0, 

°y  m  -  *8^  £(<  "  1)  Px  +  <*  +  3>  Py]  U7> 

so  o^  *  -  p^.  Equations  (16)  and  (17)  do  not  give  a  normal  stress  equal 

to  the  far-field  stress  to  be  measured.  The  normal  stress  is  Influenced 

by  the  Poisson's  ratio  for  the  medium  and  the  loading  ratio  p_/pv*  From 

*  y 


Equation  (17),  estimates  of  the  measurement  error  can  be  obtained  if  we 
examine  the  two  numerical  examples: 


v  •  1/4  (ie  »  2) 

v  -  1/3  (k  -  5/3)  ' 


If  p_/pv  <  1,  the  errors  are  acceptable. 

Although  the  gage  compliance  is  omitted  from  the  Initial 
considerations,  the  indications  are  that  slip  on  the  gage  faces  should 
be  encouraged. 


67 


SECTION  5 


NUMERICAL  SIMULATION  OF  THE  DYNAMIC  RESPONSE  OF  A  FOIL  INCLUSION 

5.1  INTRODUCTION 

He  have  numerically  simulated  the  response  of  a  thin  foil  ytterbium 
gage  in  a  PMMA  matrix  subjected  to  a  shock  wave.  Figure  26  shows  the 
configuration  of  experiments  performed  by  Y.  M.  Gupta. ^  The  matrix  is 
composed  of  four  parts:  a  1-mm  cover  plate,  two  side  pieces,  and  a 
bottom  piece.  The  shock  wave  is  induced  by  a  3-mm-thlck  flyer  plate 
that  impacts  the  cover  plate.  Four  foils  are  embedded  in  the  PMMA:  two 
with  the  foil  plane  parallel  to  the  shock  front  and  two  with  planes 
perpendicular  to  the  shock  front.  The  present  discussion  concerns  one 
of  the  foils  that  is  parallel  to  the  shock  front.  Figure  27  shows 
details  of  the  ytterbium  foil. 

A  peak  shock  wave  stress  in  the  PMMA  of  8.5  kbar  is  induced  by  the 
flyer  plate.  The  tool  for  modeling  the  gage-matrix  interaction  is  a 
two-dimensional  code;  thus,  we  must  reduce  the  problem  to  a  plane 
symmetry  case.  This  is  done  by  assuming  that  the  distance  between  the 
legs  of  the  foil  is  large  compared  with  the  foil  width  and  the  pulse 
length.  Figure  28  shows  the  area  modeled  in  the  computation.  The 
interface  between  the  foil  and  the  matrix  is  a  slide  line  with  no 
friction.  Separation  of  the  foil  and  matrix  is  allowed  to  occur.  The 
left  boundary  is  at  the  center  of  the  foil  width  and  forms  a  symmetry 
boundary.  The  right  boundary  of  the  matrix  will  be  a  symmetry  boundary 
until  reflected  signals  from  the  gage  reach  the  boundary.  The  right 
boundary  is  modeled  as  a  symmetry  boundary  all  through  the  calculation. 
This  right  boundary  is  considered  far  enough  from  the  foil  for  the  foil 
stresses  to  be  negligibly  affected  by  the  assumption  of  symmetry.  The 
lower  boundary  is  modeled  as  a  continuum  boundary.  The  mathematical 
model  in  the  code  for  a  continuum  boundary  is  not  exact,  but  it  gives  a 
small  reflection  for  a  plane  shock  wave  at  normal  incidence.  Note  that 
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the  signals  fro*  both  ths  right  and  lower  boundaries  arrive  at  the  gage 
later  than  the  rarefaction  from  the  top,  and  therefore  do  not  Influence 
the  foil  response.  The  Interface  between  the  target  and  the  flyer  plate 
la  a  slide  line  with  no  friction.  Finally,  the  top  flyer-plate  boundary 
Is  a  free  surface  and  Is  modeled  as  such. 

Figure  29  shows  the  mesh  In  the  neighborhood  of  the  gage.  The 
difference  In  cell  slse  between  the  ytterbium  and  the  PMMA  should  give  a 
correct  Impedance  match  of  the  materials. 

5.2  MATERIAL  MODELS 

The  constitutive  model  for  FNMA  Is 

P  -  Kj  H  +  Kj  |i2  +  (y  -  1)E  p  (18a) 


where 


P 

P 

P 

*1 

*2 

E 


and 


’ll 


where 


pressure 

p/Po  -  1 

density  (PQ  m  Initial  density) 
0.07243  Mbar 
0.17515  Mbar 
Internal  energy  density 
2.0 


2  G  e 


1J 


stress  deviator  rate 
strain  deviator  rate 


(18b) 


Shear  modulus  G  ■  0.025  Mbar 

Von  Mises  yield  model  with  TQ  ■  0.0045  Mbar 

Initial  density  p  “  1.184  g/cm^. 


R  (cm  x  10“3) 


JA-6877-9 


Figure  29.  Mesh  in  the  neighborhood  of  the  foil. 


Soae  uncertainty  exist*  about  tha  proper  yield  value  of  PMMA.  It  seens 
reasonable  to  believe  that  the  present  value  is  soaewhat  too  high. 

Gupta  has  obtained  the  stress  deviator  in  the  shock  wave  direction 

O 

versus  pressure  for  a  uniaxial  strain  experlaent.  The  stress  deviator 
peaks  at  2  kbar  for  a  pressure  of  7  khar  and  then  decreases  because  of 
strain  softening*  If  we  assuae  that  the  aaterial  yields  at  that  point, 
a  slaple  calculation  shows  that  the  yield  value  is  3  kbar. 

The  constitutive  aodel  for  ytterbiua  obtained  from  Gupta  has  a 
different  fora  than  the  one  for  PMMA.  The  PMMA.  yield  aodel  refers  the 
yield  value  to  the  yield  in  slaple  tension  as  opposed  to  the  aodel  for 
ytterbiua,  which  refers  yield  to  the  second  Invariant*  The  conversion 
factor  is  3°**.  The  ytterbiua  aodel  was  converted  to  the  fora  for  PMMA 
using  y  “  2.0,  Kj  ■  0.148  Mbar,  ^  •  0.0,  and  an  initial  density  of  6.98 
g/ca^.  The  stress-strain  deviator  relationship  is  siailar  with  a  shear 
aodulus  of  0.0725  Mbar,  but  the  flow  rule  is  different: 


(ip0*5  -  Tq  -  MEP  <  0 


where 

H 

- 

°.5  .tJ  .tJ 

K* 

- 

,  p  p  .0.5 
<°.5  ejj) 

eiJ 

- 

coaponent  of  the  plastic  strain 

Y 

a 

0.000473  Mbar 

0 

M 

• 

0.00454  Mbar. 

Q 

The  plezoreslstant  aodel  for  ytterbiua  is  given  by  Gupta.  After 
specialising  to  the  present  syaaetry,  we  have 


&s  -  1U  +  Ra  (A3P  +  B2«  +  rfi p) 

0  0  zz 


(19) 


where 

**o  " 


inital  resistivity 


P  ■  pressure 

a  ■  stress  In  the  direction  normal  to  the  plane  of 
symmetry 

EP  ■  square  root  of  the  second  invariant  of  plastic  strain 

A  •  7.84  per  Mbar 

B  "  16.43  per  Mbar 

17  ■  1.13  diaensionle8S. 

The  resistance  (per  unit  length)  for  a  single  cell  is  R  ■  Rs/Ac,  where 
Ac  “  current  area  of  cell.  Similarly,  the  initial  resistance  is  RQ  - 
Rs0/A0.  We  note  that  EP  in  Equation  (19)  represents  contributions  from 
the  Increase  in  defects  and  the  like;  dimensional  changes  are  accounted 
for  by  the  area  term  in  the  formula  for  the  resistance.  For  the  plane 
symmetry  assumed  here,  Kirchoff's  formula  for  the  summation  of  parallel 
resistance  is  valid  (see  Appendix):  1/R  “  1/R^  +  I/R2  ...  +  1/RQ. 
Assuming  RsQ  to  be  the  same  for  all  the  cells,  we  obtain  after  some 
algebraic  manipulations,  the  following  formula  for  the  relative  resis¬ 
tance  change  due  to  the  total  cross  section  change  of  the  ytterbium 
foil: 

(R  -  Rq)/Ro  -  Z  Ao/{Z  (Ac/Fr)}  -  1  (20) 

Where  E  stands  for  summation  over  the  cells,  and  Fr  “  Rs/Rs0.  It 
is  unnecessary  to  express  the  dimensional  changes  with  strains  because 
the  cell  area  is  directly  available  in  the  code  calculation.  Also,  it 
is  not  necessary  to  know  the  Initial  resistivity. 

5.3  RESULTS 

The  simulation  was  performed  to  3.38  ps.  The  rarefaction  from  the 
top  of  the  flyer  plate  arrives  at  about  2  pa.  The  stresses  have  in  no 
way  equilibrated  (l.e.,  generally  gone  to  zero)  when  the  calculation 
stops,  but  they  have  decreased  nearly  two  orders  of  magnitude  compared 
with  the  initial  shock  wave  stress  of  8.5  kbar.  Figure  30  shows  the 
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material  boundaries  and  mesh  around  the  foil  at  3.38  ps.  The  gage  has 
separated  from  the  matrix  material  because  the  inclusion  area  in  the 
matrix  is  larger  than  it  was  initially.  To  see  why  this  happens,  it  is 
necessary  to  consider  the  motion  when  the  initial  shock  wave  hits  the 
gage.  The  ytterbium  is  about  six  times  denser  than  the  PMMA.  This 
density  ratio  means  that  the  particle  motion  of  the  matrix  material 
above  the  gage  decreases,  producing  a  bend  in  the  matrix.  No  separation 
appears  at  this  time  because  the  material  below  the  gage  and  the  gage 
Itself  expand  to  fill  the  void  as  a  result  of  overpressures.  If  the 
yield  strength  of  the  matrix  is  low  enough,  the  material  above  the  gage 
and  immediately  to  the  right  will  deform  plastically,  and  the  deforma¬ 
tion  will  be  frozen  in  when  the  stresses  are  removed.  This  qualitative 
analysis  appears  reasonable,  considering  that  the  yield  stress  of  PMMA 
is  4.5  kbar  and  the  shock  wave  stress  Is  8.5  kbar. 

To  study  separation  and  shear  at  the  interface  between  the  foil  and 
the  matrix,  we  followed  the  motion  of  Interface  nodes  that  were  at  the 
center  part  of  the  foil.  The  nodes  were  initially  in  the  same  position: 
one  pair  on  the  upper  Interface  and  one  pair  on  the  lower  Interface. 

The  results  indicate  that  no  significant  separation  or  shear  occurred 
until  the  stresses  started  to  relax. 

We  have  sampled  Lagrangian  time  histories  for  the  whole  ytterbium 
foil  and  three  single  cells  (marked  with  an  asterisk  in  Figure  28)  in 
the  foil.  For  the  whole  foil,  all  quantities  except  the  relative 
resistance  change  are  averaged  without  weighing.  The  relative 
resistance  change  is  calculated  with  Equation  (20).  Figure  31  shows 
relative  resistance  change  versus  time  for  the  whole  foil.  The  peak 
value  la  68Z,  and  the  residual  at  3.38  }is  is  14Z.  The  experimental 
values  (from  uniaxial  strain  data)  are  75Z  and  15Z,  respectively. 

Figure  32  shows  stress  in  the  z  direction  for  the  whole  foil  versus 
time.  From  the  three  single  cell  histories,  it  is  clear  that  the 
residual  resistance  varies  significantly  between  the  cells,  whereas  the 
z  stress  is  almost  zero,  as  for  the  whole  foil.  The  resistance  curve  in 
Figure  31  is  therefore  not  typical  for  the  individual  parts  of  the  foil. 
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In  particular,  the  right  end  of  the  foil  has  a  wildly  varying  resistance 
This  is  due  to  the  large  amount  of  plastic  deformation  caused  by  the 
complex  stress  state  there. 

5.4  DISCUSSION 

There  are  several  reasons  for  the  differences  in  resistance  change 
at  peak  and  unloading  between  computation  and  experiment.  Possible 
sources  of  errors  are  constitutive  models  (for  ytterbium  and  PMMA)  and 
the  resistivity  model  of  Equation  (19);  that  is,  the  constants  in  the 
models  or  the  models  themselves  may  be  in  error.  A  significant  differ¬ 
ence  appears  to  be  due  to  the  plastic  strain  term  in  Equation  (19). 
Figure  33  shows  average  Ep  for  the  foil  versus  time.  The  final  value  is 
0.24;  the  relative  resistance  change  is  14%.  With  q  ■  1.13  in  Equation 
(19),  one  would  expect  a  residual  of  about  27%.  This  disagreement  is 
explained  by  a  residual  change  in  the  cross  section  area. 

Figure  34  presents  average  density  for  the  foil  versus  time.  As 
shown,  a  net  expansion  of  about  7%  is  obtained,  which  for  our  plane 
symmetry  implies  the  same  Increase  in  area.  A  certain  residual  stress 
also  exists,  but  the  main  effect  comes  from  the  residual  in  the  area. 
Using  the  resistivity  formula  (19),  we  can  calculate  the  relative 
resistivity  Ff  -  \/*BO  "l  +  n*Ep-l  +  1.13  *  0.24  »  1.27,  where  the 
stresses  are  assumed  to  be  zero.  The  relative  resistance  change  for  the 
oil  is  from  (20)  (R  -  R0)/RQ  -  V<Ac/Fr>  -  1  -  1/1.07  *  1.27  -  1  -  19%. 
This  does  not  agree  with  the  above  value  of  14%  that  was  calculated  by 
using  electrical  circuit  theorems  viewing  each  cell  as  an  electrical 
wire;  however,  this  lack  of  agreement  is  reasonable  because  the  use  of 
average  values  is  simply  not  accurate.  It  must  be  noted  that  the 
plastic  strain,  and  thereby  the  resistivity,  varied  considerably  along 
the  cross  section  of  the  foil. 

More  troublesome  is  the  comparatively  large  increase  in  cross 
section  area.  From  the  pressure  relationship  [Equation  (18a)),  it  is 
clear  that  a  7%  Increase  in  area  would  generate  a  tension  of  several 
kilobars,  which  is  not  found  in  the  actual  calculation.  Originally,  it 


was  assumed  that  this  could  be  explained  by  the  Internal  energy  tern  in 
Equation  (18a)— that  is,  the  Internal  energy  increases  as  a  result  o£ 
the  plastic  work;  however,  a  hand  calculation  shows  that  the  Internal 
energy  is  at  least  one  order  of  magnitude  too  snail*  The  reason  for 
this  discrepancy  is  not  clear.  One  possible  explanation  is  that  the 
Increase  In  area  cannot  be  reliably  estimated  by  calculating  the  average 
density.  A  similar  point  nay  then  hold  for  our  calculation  of  0.24  for 
the  average  plastic  strain.  Future  calculations  will  address  this 
problem. 

The  plastic  strain,  as  currently  defined,  does  not  necessarily 
increase  monotonically.  The  average  plastic  strain  in  Figure  30  is 
nonotonlc,  but  this  is  not  true  for  a  single  cell.  The  plastic  strain 
component  in  Equation  (19)  includes  a  contribution  from  defect  produce 
tion  in  the  material.  One  expects  the  number  of  defects  to  Increase 
monotonically  independent  of  the  direction  of  the  deformation.  Both 
plastic  work  and  Hill's  definition  of  equivalent  plastic  strain  are 
strictly  monotonic.  After  completing  our  calculations,  we  found  that 
Gupta  does  indeed  use  Hill's  definition  for  deriving  the  plastic  strain 
coefficient  q  in  (19).  We  will  use  Hill's  definition  of  equivalent 
plastic  strain  in  future  calculations. 


5.5  CONCLUSIONS 

The  calculation  showed  reasonable  agreement  with  experiment.  The 
calculated  relative  resistance  change  peaked  at  68Z  versus  75Z  for  the 
experimentally  observed  change.  The  residual  was  14Z  by  calculating  and 
15Z  by  experimental  observation. 

The  residual  resistance  Increase  depends  not  only  on  the  plastic 
term  in  Equation  (19),  but  also  on  a  permanent  Increase  in  the  area  of 
cross  section  of  the  foil. 
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SECTION  6 


EXPERIMENTAL  INVESTIGATIONS 

6.1  INTRODUCTION 

The  experimental  Investigations  address  the  effects  on  stress 
measurement  of  the  inclusions  formed  by  the  stress  gage  package  and  by 
the  borehole  fill.  Two  types  of  experiments  are  envisaged:  (1)  static 
trlaxlal  loading  of  a  block  of  geologic  simulant,  with  and  without  a 
borehole,  and  (2)  dynamic  uniaxial  strain  loading  of  matrices  containing 
gage  packages.  They  will  provide  known  free-field  stresses  for 
comparison  with  gage  measurements.  They  are  also  designed  to  provide 
results  for  comparison  with  theoretical  predictions  so  that  the 
mechanics  of  stress  measurement  can  be  studied. 

We  have  designed,  constructed,  and  proof-tested  the  apparatus  for 
the  static  experiments.  Rather  than  proceed  with  experiments  to  show 
the  effect  of  borehole  fill,  we  have  examined  sensors  suitable  for 
providing  experimental  results,  such  as  slip  at  borehole  flll/medla 
Interfaces.  This  work  is  discussed  in  Section  6.3. 

The  design  of  the  dynamic  experiment  is  the  result  of  calculations 
to  determine  the  feasibility  of  using  a  4-inch-diameter  gas  gun  and 
reduced-scale  flatpack  gages  to  investigate  the  relationship  between 
free-field  and  gage  stresses.  We  found  that  gage-medium  equilibrium 
could  not  be  achieved  during  pulse  durations  obtainable  from  4-inch- 
diameter  gas  gun  flyer  plates  unless  the  steel  armor  of  the  gage  were 
reduced  to  very  thin  foils.  Because  this  reduction  would  alter  gage 
response,  we  concluded  that  longer  pulses  than  are  achievable  in  the 
4-lnch-dlameter  gas  gun  would  be  required.  An  experiment  that  fulfills 
the  requirements  is  described  in  Section  6.4. 


6.2  STATIC  EXPERIMENTS 

Figure  35  shows  the  experimental  apparatus  designed  to  apply  a 
general  trlaxial  stress  field  to  a  rectangular  block  of  geologic 
material  or  rock  slaulant.  The  block  (10  x  10  z  22  Inches)  is  confined 
by  Inserts  within  a  pressure  vessel.  The  pressure  Is  applied  by  three 
pairs  of  flstjacks  on  the  faces  of  the  block,  each  pair  applying  the 
sane  pressure*  As  shown,  the  block  can  have  a  central  borehole  fill 
Inclusion  (2  Inches  In  diameter).  The  apparatus  Is  designed  for  a 
maximum  pressure  difference  of  0.5  kbar  and  a  hydrostatic  pressure  of 
1  kbar.  It  can  also  provide  plane  strain  loading  for  correlation  of 
results  with  the  theoretical  predictions  of  the  finite-element 
technique.  The  hydraulic  system  was  tested  to  2000  psi  on  a  block  of 
rock  simulant.  Flaws  in  the  flatjack  manufacture  were  corrected. 

Three  basic  experiments  are  planned  with  rock  simulant  blocks: 

(1)  no  borehole,  (2)  borehole  with  unbonded  grout  fill,  and  (3)  borehole 
with  bonded  grout  fill.  We  are  currently  experimenting  with  miniature 
stress  and  strain  gages  suitable  for  embedding  In  grout  while  minimizing 
the  Inclusion  effects  of  the  gages  themselves.  Sensing  elements  In  thin 
fluid  packages  appear  to  be  suitable  candidate  gages. 

6.3  GAGE  DEVELOPMENT  FOR  LABORATORY  EXPERIMENTS 

Two  types  of  gages  are  being  developed  for  use  In  our  static 
borehole  Inclusion  experiments.  The  first  type  Is  a  miniature  stress 
gage  designed  to  minimize  the  Inclusion  effect  on  stress  measurement  In 
small-scale  laboratory  experiments.  The  gage  is  much  smaller  than  the 
smallest  practical  flatpack  stress  gage.  It  also  has  a  strong  potential 
for  measuring  stress  In  dynamic  laboratory  experiments.  The  second  type 
of  gage  Is  a  small  Inductive  or  capacitive  sensor  to  measure  the 
separation  and  relative  slip  at  an  Interface,  such  as  the  Interface 
formed  by  borehole  fill  and  the  surrounding  medium.  This  gage  has  the 
potential  to  be  modified  for  dynamic  measurements. 

! 
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Hydraulic  Flatjacks. 


4  Steel  Rings 

^4  Aluminum  Inserts 

Grout  Specimen 
V/  <10  in.  x  10  in.  x  22  in.) 


(a)  Section  A-A 


Hydraulic  Lines 


(b)  Vertical  Cross  Section 


Figure  35.  Inclusion  apparatus. 
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Miniature  Stress  Gage 

The  sensing  element  of  the  stress  gage  Is  a  0.3-mil-thick,  1/4-inch 
by  1/4-inch  grid  of  ytterbium  bonded  to  a  2-mil-thick  layer  of  Kapton* 
The  leads  are  1-mil-thick  strips  of  copper  also  bonded  to  the  Kapton. 
These  elements  were  specially  prepared  for  us  by  Dynasen.*  The  resis¬ 
tance  is  about  50  ohms.  The  element  is  located  at  the  aidplane  of  a 
thin  disk  of  fluid  or  fluid-like  material.  The  disk  is  thin 'enough  to 
be  sensitive  only  to  the  stress  perpendicular  to  the  disk  face. 

C  rrently,  the  overall  dimensions  of  the  disk  are  1  inch  for  the 
diameter  and  1/10  inch  for  the  thickness,  but  these  can  be  reduced  once 
a  satisfactory  design  has  been  achieved.  The  fluid  pressure  is  approx¬ 
imately  equal  to  the  stress  to  be  measured,  and  the  change  in  resistance 
of  the  ytterbium  sensor  determines  the  fluid  pressure.  For  our 
ytterbium,  the  hydrostatic  calibration  was  6Z/kbar  for  the  resistance 
change. 

Two  kinds  of  stress  gages  have  been  tested,  the  main  difference 
between  the  designs  being  the  fluid.  In  one  design,  we  used  glycerol 
enclosed  by  urethane  sheets,  whereas  in  the  other  we  used  the  fluid-like 
material,  GE  silicone  (RTV).  Results  and  observations  from  the  testing 
to  date  of  these  two  gages  are  summarized  as  follows: 

•  Under  hydrostatic  pressure  the  gages  in  water  follow  the 
hydrostatic  calibration  curve  for  the  ytterbium  element. 

e  The  RTF  gage  measured  the  axial  stress  well  in  alluvium 
subjected  to  static  uniaxial  strain  loading,  unlike  the 
glycerol  gage. 

•  Both  types  of  gages  measured  stress  changes  satisfactorily 
above  1500  pal  (10.33  MPa)  In  2C4  grout  subjected  to 
static  uniaxial  stress  loading.  Below  1000  psl  (6.89 
MPa),  both  types  of  gage  provided  measurements  that  were 
unacceptably  low. 

•  The  RTV  gage  is  much  simpler  to  manufacture  and  is 
adaptable  for  installation  in  geologic  specimens. 


*Dynasan,  Inc.,  20  Dean  Arnold  Place,  Goleta,  California  93017. 
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In  conclusion,  the  RTv  version  of  the  stress  gage  was  better  than 
the  glycerol  version  in  performance,  ease  of  fabrication,  and  adapta¬ 
bility  to  geologic  materials.  The  unsolved  problem  that  remains  is  the 
measurement  of  stress  below  about  1500  psi.  It  is  desirable  to  have 
this  capability  so  that  the  gages  can  be  incorporated  in  the  static 
plane-strain  borehole  experiment.  If  this  problem  is  not  solved  within 
a  reasonable  time,  we  will  measure  strain  only  us^ng  miniature  rod 
gages*  embedded  and  anchored  in  the  2C4  grout. 

Because  the  RTV  version  of  the  stress  gage  showed  more  promise  than 
the  glycerol  version,  we  continued  development  work  only  on  the  RTV 
version  to  address  two  problems.  The  first  problem  occurred  when  the 
sand  or  grout  loading  caused  distortion  of  the  Inclusion  that  Involved 
stretching  of  the  midplane.  Although  the  RTV  has  a  low  shear  modulus, 
we  found  that  the  combination  of  high  pressure  and  distortion  caused 
stretching  of  the  ytterbium  grid  by  shear  traction.  The  Induced  tension 
dominated  the  resistance  change  (a  decrease).  To  suppress  this  induced 
tension,  we  bonded  the  Xapton  tape  under  the  ytterbium  grid  to  a  thin 
steel  plate.  The  resulting  gage  response  was  a  correct  measurement  of 
stress  change,  except  again  for  an  insensitivity  at  low  stresses  «  1000 
psi). 

The  second  problem  is  the  applicability  of  the  plezoresistive 
constants  that  have  been  determined  for  ytterbium  foil,  to  ytterbium 
that  has  been  vapor  deposited,  as  is  the  case  with  our  grids.  We  solved 
this  problem  by  hydrostatically  loading  two  RTV  gages  to  1  kbar.  In  one 
gage,  the  ytterbium  grid  was  free  and  in  the  other  the  ytterbium  grid 
was  bonded  to  a  thin  steel  plate.  The  resistance  change  was  calculated 
using  the  plezoresistive  coefficients  for  ytterbium  foil  and  the  results 
were  compared  with  the  experimental  values.  The  extent  of  the  agreement 
is  shown  in  Figure  36,  the  calculated  values  of  the  linear  pressure- 
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resistance  change  curves  being  6.0Z/kbar  and  4.7Z/kbar  for  the  free  and 
bonded  grids.  We  conclude  that  the  piezoresistive  coefficients  are 
applicable  to  our  grids  in  the  pressure  range  of  our  investigation  (<  1 
kbar).  We  have  yet  to  resolve  the  observation  that  the  loading  and 
unloading  lines  in  Figure  36  are  not  the  same. 

Interface  Gage 

Interface  gages  are  based  on  the -principle  that  inductive  loops  or 
capacitive  plates  located  on  either  side  of  an  interface  will  change 
their  mutual  inductance,  M,  or  capacitance,  C,  when  the  interface 
experiences  separation  or  relative  slip.  These  gages  form  part  of  a 
tank  circuit  that  oscillates  at  a  design  frequency.  When  M  or  C 
changes,  the  natural  frequency  changes.  Design  calculations  Indicate 
that  the  mutual  inductance  gage  is  more  sensitive  that  the  capacitance 
gage  because  compact  ferritic  cores  with  many  turns  of  wire  are  avail¬ 
able  to  produce  high  mutual  Inductances.  The  resulting  gage  is  about 
1/4  inch  in  diameter  and  1/8  inch  thick.  The  high  Inductance  allows  us 
to  use  standard  electronic  equipment  and  frequencies  below  1  MHz. 


6.4  LONG-DURATION  PLANE  WAVE  EXPERIMENT 

One  objective  of  the  dynamic  experimental  program  is  to  perform 
well-controlled  uniaxial  strain  Impact  experiments  on  matrix  materials 
containing  embedded  flatpack  gages  with  plezoresistant  foil  sensors. 

The  results  of  these  experiments  are  Intended  to  be  compared  to  those 
from  analogous  gas  gun  experiments  performed  by  Gupta,  in  which  the 
matrix  contained  bare  plezoresistant  foils  without  the  surrounding  gage 
package. 

The  SRI  4-lnch  gas  gun  would  be  a  convenient  loader  for  these 
experiments,  if  the  valid  one-dimensional  test  time  (before  the  arrival 
of  edge  effects)  is  great  enough  for  a  scaled  down  flatpack  gage  package 
to  come  to  mechanical  equilibrium  with  the  flow  in  the  surrounding 
matrix  material.  For  a  PMMA  matrix  material,  the  one-dimensional  test 
time  in  the  4-lnch  gas  gun  is  estimated  to  be  about  8-9  ps  over  the 


central  Inch  of  the  target  at  the  stresses  of  Interest.  Therefore,  we 
performed  a  series  of  calculations  with  the  SRI  PUFF  one-dimensional 
finite  difference  code  to  determine  the  response  time  of  thin  gage 
packages  to  uniaxial  strain  shock  compressions  and  release. 

For  these  one-dimensional  calculations,  the  gage  package  is  modeled 
simply  as  two  steel  layers  surrounding  a  PMMA  layer  (at  the  center  of 
which  would  be  'the  plezoreslstant  foil).  The  package  is  embedded  in  a 
PMMA.  matrix,  and  the  matrix  is  impacted  by  a  PMMA  flier.  The  problem  is 
run  for  8  ps,  the  estimated  maximum  valid  one-dimensional  test  time  in 
the  4- inch  gas  gun.  tfe  performed  the  following  series  of  calculations. 

First,  we  performed  thick  (1-inch)  flier  calculations  to  determine 
the  response  to  compression  only  (no  rarefaction  from  the  free  surface 
of  the  impactor  for  at  least  8  ps)  for  two  gage  package  thicknesses. 

The  equilibration  time  at  the  center  of  the  package  to  the  compression 
was  determined  and  the  appropriate  flier  thickness  was  calculated.  We 
then  performed  thin  flier  calculations  to  determine  the  response  time  of 
these  gage  packages  to  step  compression  and  release.  We  considered 
three  flier  velocities,  0.56,  1.0,  and  1.56  km/s  (the  latter  being  the 
maximum  velocity  achievable  in  this  gun),  corresponding  to  PMMA  stresses 
of  about  1,  2,  and  about  3.7  GPa. 

We  verified  that  the  computational  results  were  not  zone  size 
limited  by  halving  the  zone  size  in  one  calculation.  We  established 
baseline  computational  rise  times  and  numerical  oscillation  charac¬ 
teristics  for  use  in  Interpreting  the  gage  package  calculations  by 
performing  calculations  with  the  same  zoning  but  with  all  components 
(target,  gage  package,  and  impactor)  PMMA. 

We  computationally  examined  two  gage  package  thicknesses.  In  the 
thicker  package,  each  steel  layer  was  0.02  inch  thick  and  the  PMMA  was 
0.10  inch  thick.  In  the  thinner  package,  each  of  these  dimensions  was 
halved.  Although  even  thinner  packages  could  conceivably  be  fabricated, 
the  steel  layers  would  be  foils  no  longer  having  the  characteristics  of 
the  full-scale  components  and  the  results  would  not  be  relevant  to  real 
gage  packages. 


The  calculations  led  to  two  primary  conclusions*  First,  only  the 
congressional  response,  not  the  following  release,  could  be  examined  on 
the  4-inch  gas  gun.  Because  the  unloading  response  is  equally 
important,  the  limited  duration  of  planar  flow  in  the  4-inch  gas  gun  is 
a  severe  restriction.  Second,  even  in  compression,  only  the  thinner 
package  reaches  mechanical  equilibrium  rapidly  enough  to  be  tested  in 
the /4-inch  gun. 

These  conclusions  pertain  only  to  equilibration  in  the  wave 
propagation  direction  (one-dimensional  equilibration).  The  actual 
package  would  have  to  equilibrate  in  the  transverse  directions  as 
well.  Although  the  equilibration  times  for  these  processes  are  longer, 
the  magnitude  of  the  nonequilibrium  effects  are  smaller. 

Suitable  stress  pulses  for  measuring  both  loading  and  unloading 
response  in  realistic  embedded  gage  packages  can  be  generated  by  an 
explosive  launcher  device  that  propels  flyer  plates  much  larger  in 
diameter  and  thickness  than  those  used  with  gas  guns.  The  design  of 
such  a  launcher  is  presented  below.  It  produces  plane  wave  pulse 
durations  long  enough  for  an  examination  of  the  loading,  peak  stress, 
and  unloading  responses  of  stress  gage  packages. 

The  design  presented  is  specific.  Variations  to  extend  the  stress 
or  duration  range  are  summarized  immediately  following  the  design. 


Flyer  Plate/Target  Characteristics 

We  choose  the  plastic  PMMA  as  the  flyer  plate  material  because  it 
is  light  and  the  properties  are  known.  We  choose  the  rock  simulant  2C4 
for  the  target  material  because  the  stress  gage  package  can  be  cast  in 
place  and  again  the  properties  are  known.  For  planar  impact  at  a 
velocity  small  compared  with  shock  velocities,  the  stress  Induced  in  the 
target  is  approximately 


a 


plcl  p2c2 
plcl  +  P2C2 


(21) 
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fear's* 


!*HP 


where 

p^  ■  flyer  plate  density 

c^  ■  flyer  plate  shock  velocity 

p2  *  target  density 

C2  ■  target  shock  velocity 

V  “  Impact  velocity 

a  ”  stress  In  target. 

For  a  flyer  plate  of  PMMA  and  a  target  of  2C4,  we  have 


3  3 

Pj  ■  1.18  g/cm  Cj^  ■  3.5  mm/ps  P2  "  2.15  g/cm  c2  ”  3.3  nm/ps 


and  when  V  *  100  m/s,  the  stress  In  the  target  Is 


a  ■' 2.6  kbar. 


(22) 


> 

& 


The  particle  velocity  In  the  2C4  Is 
plcl 


plcl  +  p?c 


V  ■  0.37V  *  37  m/s. 


(23) 


2  2 


For  a  flyer  plate  of  PMMA,  the  stress  pulse  duration,  T,  Is  twice 
the  transit  time  through  the  thickness,  h^;  that  Is,  T  »  2hj/cj.  If 
hj  -  4  in.  (10.16  cm). 


58  ps. 


(24) 


The  minimum  diameters  of  flyer  plate  and  target  are  determined  by 
the  duration  of  the  stress  pulse.  To  allow  examination  of  the  stress 
gage  package  response  to  known  unloading,  rarefaction  waves  from  the 
edge  of  the  flyer  and  target  plates  must  not  reach  the  gage  before  the 
pulse  ends.  If  the  gage  Is  located  on  the  common  axis  of  the  plates 
near  the  Impact  surface  and  if  the  velocity  of  the  rarefaction  wave  is 
approximated  by  the  shock  velocity,  the  minimum  diameter  Is  D  ■ 

2  max(ci,  For  a  4-inch-thick  PMMA  flyer  plate  striking  a  2C4 
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target,  we  hare  max  (c^,  c^)  ■  ■  3.5  mm/ps  and  T  »  58  ps,  so  min  D  » 

16  In.  (40.6  cm) • 

In  our  design  we  choose  D  -  21  inches  to  ensure  that  the  arrival 
ties  of  the  edge  rarefaction  (l.e.,  76  (is)  is  well  past  the  pulse 
unloading  tlae. 

The  minimum  thickness  of  the  2C4  is  similarly  chosen  to  allow  an 
adequate  unloading.  With  5  inches  of  2C4  behind  the  gage  package,  the 
rarefaction  from  the  near  face  of  the  target  will  arrive  at  77  pa. 


Launcher  Dynamics 

Figure  37  shows  schematically  the  launcher,  flyer  plate,  and 
target.  A  steel  pressure  vessel  is  filled  with  a  stoichiometric  mixture 
of  hydrogen  and  oxygen  at  some  specified  pressure,  and  the  mixture  is 
held  by  a  diaphragm.  Detonation  of  this  gaseous  mixture  Increases  the 
pressure,  which  causes  the  diaphragm  to  rupture,  and  the  expanding  steam 
drives  a  piston  with  an  attached  flyer  plate  along  a  short  tube  toward  a 
target. 

The  pressure  vessel  with  its  diaphragm  arrangement  is  the  vessel 
used  to  launch  turbine  missiles  at  reinforced  concrete  barriers  in  our 
past  industrial  safety  experiments.1^  We  found  this  energy  source  to  be 
reliable,  reproducible,  versatile,  and  economical. 

Our  design  is  based  on  launching  a  16-inch-diameter  piston  capable 
of  carrying  a  21-inch-diameter  flyer  plate  of  4-inch-thick  PMMA.  The 
flyer  plate  nay  be  thicker  and  larger  in  diameter  if  longer  pulses  are 
required* 

Let  the  piston  assembly  consist  of  three  main  parts:  a  16-inch- 
dlameter,  2-inch-thick  aluminum  pusher  plate;  a  16-lnch-dlaneter, 
1/2-inch-thick,  20-lnch-long  aluminum  tube;  and  the  21-lnch-diameter, 
4-inch-thick  PMMA  flyer  plate.  The  masses  of  these  parts  are  approx¬ 
imately  17.8,  22.2,  and  26.8  kg,  respectively,  giving  a  total  piston 
mass  of  66.8  kg.  We  use  a  mass  of  M  ■  70  kg  in  our  design  calculations 
to  take  into  account  any  additional  mass  from  attachments,  overslslng, 
and  the  like. 


If  we  design  for  an  impact  velocity  of  V  ■  100  m/s,  the  piston 
kinetic  energy  to  be  supplied  by  gas  work  is 


12  5 

K  -  j  Wf  -  3.5  *  10 


for  M  -  70  kg.  An  aquation  of  state  appropriate  for  steam  at  the 
pressures  in  our  launcher  «  300  bar)  is  the  perfect  gas  law  pyY  - 
constant,  with  y  ■  1.14.  The  maxima -gas  work  available  when  the 
Initial  pressure  and  volume  are  p0  and  vQ  is 

w(o)  .  Vo_ 

w-  y  -  1 

Our  pressure  vessel  has  a  volume  of  vQ  -  0.0287  m^,  which,  with  a 
constant  volume  pressure  following  detonation  of  pQ  -  3250  psi  (224 
bars),  gives  the  gas  work  as 


,«»  - 


45.9  x  10  Mi 


This  work  is  more  than  13  times  the  kinetic  energy  required.  The 
pressure  of  the  hydrogen/oxygen  mixture  before  detonation  that  is 
required  to  produce  a  constant  volume  pressure  of  224  bars  is  21  bars 
(300  psi,  ratio  of  10.8).  After  the  diaphragm  splits,  the  gas  expands 
to  a  volume  of  v^  •  0.0332  m^  before  it  starts  to  move  the  piston.  The 
expansion  comprises  the  orifice  volume  (4  Inches  in  diameter,  6  Inches 
long)  and  a  1-inch  space  behind  the  piston  for  a  smoother  launch.  When 
v  -  Vj,  the  pressure  is  -  2810  psi  (194  bars),  and  the  work  done  on 
the  piston  when  the  steam  has  Increased  in  volume  from  v1  to  v2  is 

w  -  W<L)  1  -  1  (28 


where 


W(1)  .  Vi 

y  -  1 


45.9  x  10 J  Nw 


n. 


The  value  of  v2  whan  tha  gas  work,  W,  equals  the  required  kinetic  energy 
of  K  ■  3.5  x  105  h  Is  therefore 


*2  -  vt 


[l'^l 


r-i  _ 


1.76  v. 


%  3 

Because  vj  -  0.0332  r ,  we  have  V2  -  0.0584  a  ;  hence  the  Increase  In 
voluaa  caused  by  piston  notion  i*  v2  -  V1  “  0.0252  a^.  The  16-inch- 
dlaaeter  piston  cross-sectional  area  Is  A  ■  0.13  a*;  thus  the  piston 
travel  before  It  attains  a  velocity  of  100  a/s  is  L  ■  0.19  a  (7.5  Inches) 
When  the  steaa  has  expanded  to  the  voluae  v2,  the  pressure  has  decreased 
froa  pi  •  194  bars  to  p2  “  102  bars. 


An  upper  bound  estlaate  for  the  duration  of  piston  travel  Is 
obtainable  by  approxlaatlng  the  pressure  in  the  equation  of  notion 


Mx  -  pA 


P  «  Pt  -  (PL  -  P2)  l 


Integration  then  leads  to  the  expression  for  the  duration  of  piston 


notion 


_  [  ML  11/Z  -1  p2 

K  COB 


For  M  ■  70  kg,  L  •  0.19  a,  A  ■  0.13  a*,  p^  ■  194  bars,  end  p2  -  102 
bars,  we  obtain  t  ■  3.4  as. 


Structural  Pas It 


In  calculating  the  aass  .of  the  piston,  we  assuasd  that  a  2-inch- 
thick,  16-lnch-dlaaater  piston  pusher  plate  of  alualnun  alloy  would  be 
strong  enough  to  withstand  the  gas  pressure.  When  the  pusher  plate  Is 
considered  as  an  elastic,  slaply  supported  circular  plate  under  a  static 


WWW 


stress  at  the  center  is 


unlforn  pressure  p^,  the 
3(3  +  v)a2Pl 

g  m  - - - 

8h 


(32) 


where  s  is  the  radius*  h  the  thickness*  and  v  is  Poisson's  ratio.  For 
a  ■  8  Inches*  h  -  2  inches*  and  v  “  0.34,  we  have  o  -  20  p.  If  p^  -  194 
bars  (2810  psl)*  we  have  a  m  56,370  psi*  which  implies  yielding  for  most 
aluminum  alloys,  especially  since  the -load  is  applied  suddenly.  Because 
the  piston  is  not  reusable*  some  yielding  is  acceptable.  The  static 
plastic  collapse  pressure  of  the  plate  is 


3h 

2a 


2  *y 


32  °y 


(33) 


where  ay  is  the  yield  stress.  If  «Jy  -  66,000  psl  (7075-T6  alloy),  the 
static  collapse  pressure  is  6,188  psl  (427  bars),  which  is  well  above 
p|.  Thus,  without  further  analysis  we  may  say  that  a  2-inch-thick  plate 
will  probably  be  adequate.  A  3-inch-thick  plate  will  certainly  be 
adequate,  but  it  Increases  the  piston  mass  and  is  more  expensive. 

The  2 0-inch- long  tubular  portion  of  the  piston  (Figure  37)  is 
designed  to  be  longer  than  the  diameter  for  smooth  travel  along  the 
steel  launcher  tube.  The  maximum  thrust  on  the  piston,  which  occurs  at 
the  start  of  motion,  is  25.2  x  105  R.  If  the  well  thickness  is  1/2 
inch,  the  maximum  “static”  axial  stress  is  23,200  psl  (1,600  bars), 
which  is  well  below  the  yield  stress  for  most  aluminum  alloys,  and  the 
radius- to- thickness  ratio  of  16  la  small  enough  to  preclude  buckling 
under  these  loading  conditions. 

The  flyer  plate  is  stressed  by  inertial  forces,  which  are  greatest  at 
the  start  of  notion  when  the  thrust  on  the  piston  is  F  •  25.2  x  105  N. 

For  a  piston  mass  of  M  -  70  kg,  the  initial  acceleration  is  a  -  F/M  ■ 
36,000  m/s2.  A  4-inch-thick  PMMA  flyer  plate  has  an  areal  density  of 
m  ■  12  g/cn2;  hence,  the  uniformly  distributed  inertial  pressure  is  as  - 
43.2  bars  (626  psl). 
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Am  clastic  circular  plate  of  radius  r j  on  a  circular  support  ring 
of  radius  rQ  subjected  to  a  uniform  pressure  hss  Its  extremum  bending 
moments  at  ths  center  r  ■  0  and  oxer  the  support  ring  r  -  rQ.  These  are 
2  2 

•  ^ir  'i  ■,<*„>  ‘  tt4  ['i  -  <3  + »>]  (3*) 

where 

-  2(1  -  v)  +  (1  +  3v)02  -  2(1  +  v)02  An  02,  0  -  r^ 

for  our  flyer  plate  ms  have  -  10. 5  inches  and  rQ  “  7.75  inches,  so 
0  •  1.35.  For  FIMA,  Poisson's  ratio  v  ■  0.35;  thus  2.083  and 
Xj  -  (3  +  v)  ■  -1.267,  giving  ar(o)  •  7.82  pa  and  m^r^  ■  -6.76  p^. 

The  fiber  stress  a.  Is  related  to  the  bending  moment  by  cr.  ■  6  m  /h2, 

D  D  f 

where  h  Is  the  plate  thickness  (4  Inches).  Hence,  a.  (0)  -  2.93  p  and 

D  u 

a. (r  )  »- 1.79  p  .  For  an  inertial  pressure  of  p  ■  626  pel, 
boa  a 

ov(0)  -  1,834  pel  a.  (r  )  -  1,121  pel  (35) 

D  DO 

which  are  well  below  the  PMMA  tensile  strength  of  8,000  pel.  Thus,  the 
FIMA  flyer  plate  will  not  fracture  during  launching  because  of  its 
overall  structural  response. 

The  steel  tube  for  guiding  the  piston  and  containing  the  steam  is 
designed  to  withstand  an  Internal  pressure  of  pj  -  2,810  pel.  The  hoop 
stress  is  ■  a/h  p^.  If  the  tube  thickness  Is  chosen  to  be  h  -  2  Inches 
giving  a  ■  9  Inches  we  have  a  safe  working  stress  of  ■  12,645  psl. 


Overall  Piston-Target  Impact 

Assuming  ths  impact  of  the  piston  on  the  target  to  be  plastic  so 
that  they  both  move  together  as  a  single  mass  allows  us  to  estimate  the 
combined  velocity  after  impact*  For  a  6-inch-thick  2C4  target,  21 
inches  In  diameter,  the  mass  is  *  76  kg.  Our  piston  has  a  total  mass 
of  Mp  »  70  kg;  thus  the  combined  velocity  after  impact  Is 


When  the  piston  strikes  the  target  st  Vp  -  100  e/s,  we  get  Vpj  ■  49  m/s. 

Cheracteristics  of  Various  Flyer  Piste  Materials 

In  the  preceding  discussion,  flyer  plate/target  characteristics 
were  presented  for  the  specific  case  of  a  plastic  PMMA  flyer  impacting  a 
grout  target.  The  plastic  PMMA  flyer  Induced  a  stress  of  a  ■  2.6  kbar 
in  the  grout  target.  The  stress  induced  in  the  target  aay  be  extended 
by  appropriate  choice  of  flyer  plate  materials.  Table  2  compares  the 
flyer  plate/target  characteristics  of  4-inch-thick  aluminum,  grout,  and 
plastic  PMMA  flyers  impacting  a  grout  target.  Tor  an  Impact  velocity  of 
V  “  100  m/s,  the  stress  induced  in  the  tsrget  Increases  to  S  kbar  for  an 
aluminum  flyer  plate  and  3.6  kbar  for  a  grout  flyer  plate.  Note, 
however,'  that  the  aluminum  flyer  plate  produces  a  pulse  duration  that  is 
just  over  half  of  the  pulse  duration  produced  by  the  plastic  PMMA  flyer. 
A  thicker  aluminum  plate  could  be  used  to  increase  the  pulse  duration. 
Also  tabulated  are  the  launch  length  until  Impact,  L,  and  the  launch 
time  until  Impact,  T. 


Tabic  2.  Flyer  Plate/Target  Characteristic* 


4-Inch 

V 

Pulse  T 

Launch  L* 

Launch  T 

Flyer 

(kbar) 

(  •) 

(•) 

(«) 

PUMA 

2.6 

100 

58 

0.19 

3.4 

Alualnua 

5.0 

100 

32 

0.35 

5.7 

Groutb 

3.6 

100 

62 

0.29 

4.8 

*The  launch  distances  for  alualnua  and  grout  flyer  plates  can  be  reduced 
.  by  Increasing  the  Initial  hydrogen/oxygen  pressures. 

bThe  grout  flyer  uust  be  reinforced  because  of  Its  low  tensile  strength 
(500  pel). 
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Appendix 


COMPARISON  OF  RESISTANCE  AVERAGING  WITH  MECHANICAL 
AVERAGING  FOR  CALCULATING  THE  RESISTANCE  CHANGE 
OF  A  PIEZORESISTANCE  FOIL  IN  PLANE  STRAIN 

In  this  appendix,  we  examine  the  equivalence  of  calculating  the 
stress-induced  change  In  total  resistance  of  a  plezoreslstant  foil  In 
plane  strain  deformation  by  two  methods:  resistance  and  mechanical 
averaging.  In  resistance  averaging,  we  convert  the  stresses  and  strains 
Induced  In  the  foil  by  a  stress  wave  to  cell  quantities,  determine  the 
resistance  of  each  cell,  and  then  calculate  the  total  resistance  change 
from  a  network  circuit  analysis.  In  mechanical  averaging,  we  calculate 
an  average  value  over  the  foil  for  each  component  of  stress  and  strain 
and  use  these  to  calculate  the  total  resistance  change  of  the  foil. 

Both  methods  use  the  piezoresistance  formulation  of  T.  Gupta  expressing 
the  vector  resistance  In  terms  of  the  mechanical  state  variables  at  a 
point: 


where  R^  is  the  resistance  to  current  flowing  in  the  1  direction,  Rq  la 
the  Initial  value  of  this  (and  the  other)  resistance  components,  a  is 
stress,  e  is  strain,  la  a  plastic  work  function,  and  a,  0  and  q  are 
piezoresistance  coefficients. 

The  basic  problem  we  wish  to  calculate,  the  resistance  change  to 
constant  current  flow  In  the  length  direction  of  a  foil  as  It  Is  loaded 
by  a  normally  incident  planar  wave,  Is  shown  schematically  in  Figure 
A.l,  along  with  the  foil  configurations.  The  length  direction  of  the 
foil  has  been  arbitrarily  chosen  along  the  3-axls  and  the  Incident  plane 
wave  Is  propagating  In  the  2-dlrectlon.  The  figure  shows  the  two  foil 
orientations,  normal  and  transverse,  of  practical  interest  and  includes 
the  type  of  zoning  appropriate  for  plane  strain  calculations.  (For 
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notat tonal  convenience,  throughout  this  Appendix  we  have  set  the  length 
equal  to  one  unit  as  Indicated  In  the  figure.) 


By  aodellng  the  Internal  deformation  of  the  foil  as  two-dimensional 
plane  strain  (a  good  approximation,  usually,  since  real  foils  are 
normally  much  longer  than  they  are  wide  or  thick),  we  have  Imposed  the 
implicit  constraint  that  the  foil  ’‘faces'*  (cross  sections  perpendicular 
to  the  Indicated  length  dimension)  remain  parallel  to  the  1-2  plane. 

This  constraint  is  important  for  the  resistance  averaging  calculations. 
Below  we  describe  the  specific  resistance  and  mechanical  averaging 
approaches  used  for  the  special  cases  considered  in  this  appendix, 
formulate  a  particularly  simple  test  problem,  and  apply  resistance 
averaging  and  mechanical  averaging  to  the  test  problem  to  demonstrate 
and  compare  the  two  methods. 

Resistance  Averaging 

We  assume  that  (1)  steady  state  conditions  exist  for  the  current 
distributions  and  potential  distributions,  (2)  resistance  is  the  only 
source  of  potential  differences,  and  (3)  fringing  can  be  neglected. 

With  these  assumptions  and  a  plane  strain  representation,  the  "faces”  of 
the  foils  are  equlpotential  surfaces.  Therefore,  the  electrical  field 
$  and  current  density  3  vectors  are  In  the  3-direction  at  the  foil 
faces.  Since  each  "face"  Is  an  equlpotential  surface,  since  mechanical 
variations  are  not  permitted  In  the  3-direction,  and  since  the 
separation  of  the  faces  is  arbitrary,  then  the  potential  vectors  within 
each  cell  are  also  in  the  3-dlrectlon  only  and  the  cells  behave  as 
Independent  resistors  in  parallel. 

Thus,  in  resistance  averaging  the  total  resistance  in  the  3- 
dlrectlon,  BiOT,3*  of  the  conductor  is  obtained  by  first  calculating 
R^/Rq  for  each  cell,  j,  from  equation  (A.l)  using  the  calculated 
stresses  and  strains  in  the  cell  and  the  initial  resistance,  and  then 
summing  these  resistances  in  parallel  according  to 


'T°r'3 "  [ij  *»r 
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(A. 2) 


In  this  approach,  vs  determine  a  single  average  value  for  each 
strasa  and  strain  component  within  the  foil,  as  well  as  for  the  plastic 
work.  We  then  use  these  in  equation  (A. 1)  to  find  the  total  resistance 
change. 

Principal  engineering  strains  are  adequate  for  our  purposes.  We 
obtain  average  values  by  sunning  over  external  boundaries  of  exterior 
cells,  k,  as  indicated  below: 


where  X£  is  the  length  of  the  external  boundary  of  cell  k  in  the  i 
direction  and  x£  is  the  total  length  of  the  undeforaed  foil  in  the  1 
direction. 

To  obtain  average  stresses,  we  sun  external  forces  and  divide  by 
the  area  of  external  faces: 


i 


£  A! 


(A.A) 


To  obtain  the  average  plastic  work,  we  average  over  the  ness  or 
voluae  of  all  the  cells,  j: 


(A. 5) 


vhare  Rj  i*  the  ness  or  voluae  of  cell  j. 
Test  Problen 


Cons ids r  the  uniaxial  strain  dynamic  deformation  Induced  by  a 
constant  amplitude  planar  shock  in  a  foil  composed  of  just  two  cells,  as 


shown  in  Figure  A. 2.  Since  cross-sectional  dimensions  ere  arbitrary, 
this  can  be  either  a  transverse  or  normal  foil;  a  transverse  foil  is 
shown*  We  wish  to  calculate  the  resistance  by  the  two  methods  when  the 
deformational  wave  has  reached  the  boundary  between  cells. 

Without  loss  of  generality,  we  set  all  initial  linear  cell 
dimensions  and  all  initial  cell  resistances  to  unity*  Furthermore,  we 
postulate  that  the  foil  material  is  hydrodynamic  so  that  within  a  cell, 
Oj  •  o j  ■  Oy  We  also  select  it  to  have  hydroetatlc  piesoreelstlve 
response,  0  ■  rj  ■  0,  so  that  equation  (A.l)  becomes 

ji  -  [i  +  «<*!  +  *j  +  *3)]  flVTTnVo  <A-‘ 

O  12 

Resistance  Averaging  Applied  to  Test  Problem 

Consider  cell  1.  The  stresses  and  strains  are 


el  *  «3  "  °*  e2  5  € 

*1  "  *3  "  a2  3  °  ‘ 

Therefore,  by  equation  (A.6),  the  relative  resistance  in  this  cell 

is 

^  -  rJ  -  (1  +  3  a  a)  3-^7 
o 

For  cell  2; 


Therefore,  by  equation  (A. 2),  resistance  averaging  gives  the  total  foil 
resistance  in  the  3-direction  as 


t 
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*  -V 


.(1)  .(2) 


RT0T,3  “ 


(1  +  3  g  o)/(  1  -I-  e) 
f(T  +  3  a  <*)/(!  +  e)]+l 


1  +  3  a  o 

2  +  3  a  <3  +  e 


and  the  relative  resistance  as 


TOT. 3  _  1  /  1  +  3  ao 


l  t  1  +  3  80  \ 

"  7  (TTTooTT' 


(A.  7) 


Mechanical  Averaging  Applied  to  Test  Problem 


First,  consider  the  strains.  Within  each  cell  the  transverse 
strains  are  zero,  so 


ei  "  e3  "  0 


The  average  axial  strain  from  equation  (A. 3)  Is 


[X*  ♦  x\]  -  [(X*)°  ♦  (x*)°j 


C(x;)°  +  (x;)1] 


Ci  +ja  +  U 
u  + 


-  u .  +  n 


Now  consider  the  stresses.  Since  the  foil  material  Is  hydrostatic, 
the  three  principal  stresses  within  either  cell  are  equal.  However,  the 
average  stresses  over  the  whole  foil  may  not  be  equal  and  are  therefore 
treated  separately.  Consider  first  the  average  transverse  stresses, 


0j  and  0j.  By  symmetry  -  0^.  Now  evaluate  from  equation  (A. 4). 


where,  as  usual,  the  subscript  denotes  the  coordinate  axis  and  the 

1  2 

superscript  denotes  tbs  cell*  Remember  that  *  a,  ■  0, 
aJ  -  (1  +  e)  and  A^  -  1.  Therefore, 

1  +  c  - 

cl"°2T7"*3  * 


Now  consider  the  average  axial  stress,  o^*  Although  the  cells  are 
not  In  aechanlcal  equilibrium  In  this  direction,  ve  obtain,  by  foraally 
applying  equation  (A.4), 

1  .1  ^  2  .2 


1.1.  2.2 


Thus  the  average  stresses  In  the  foil  do  not  satisfy  the  hydrodynamic 
postulate  with  this  averaging  method  except  In  the  Halt  as  e  -*>0. 

Substituting  these  average  stresses  and  strains  Into  equation 
(A. 6),  we  find  that  the  total  relative  resistance  change  in  the  3- 
dlrection,  by  aechanlcal  averaging,  Is 


Vs  ■  [X  ♦  .  <5,  ♦  H  *  ;,)]  — ---3 

o  <1  +  Cj)(l  +  e2) 

-  [i  +  *  {J*  (j-y-f-D  +  f)]  YTT/t 


2  +  *  [la  i  »>  ♦  o] 


WJvW.V 


Conclusions 


Iha  simple  twrc*ll  transient  response  calculations  considered 
above  show  that  resistance  and  aeehanical  averaging  are  not  equivalent, 
in  general,  since  equations  (A. 7)  and  (A.8)  are  unequal.  Resistance 
averaging  was  derived  fro*  baaic  princlplea,  vhereae  mechanical 
averaging  was  simply  postulated  and  is  not  considered  to  be  physically 
reallatic.  Therefore,  in  plana  strain  calculation  of  the  type 
conaldered  here,  we  conclude  that  plesoresistance  effects  should  be 
calculated  by  resistance  averaging,  treating  the  current  carrying  cells 
in  the  length  direction  of  the  foil  as  independent  parallel  variable 
resistors. 
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